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Summary 


a 

The  Air  Force  has  a  need  for  better  estimates  of  a  parameter  known  as  "C 

2  n 
It  also  needs  to  know  how  it  (Cn)  varies  with  geographic  location,  season,  time  of 

day  and  weather  conditions.  This  need  arises  in  connection  with  optical  (LASER) 

2 

communication,  navigation  and  weapon  systems  applications,  Cn  is  related  to 
turbulence-induced,  index  of  refraction  variations  in  the  atmosphere. 

2 

The  main  purposes  of  this  report  are  to  explain  the  physical  nature  of  Cn,  to 
describe  the  current  state  of  the  art  in  this  area,  and  to  make  suggestions  for 
future  research. 

2 

Section  2  shows  how  the  parameter  Cn  arises  in  atmospheric  turbulence. 

Quantitatively,  Cn  =  (S')  ft  '  where  (0')z  is  the  mean  square  variation  of  the 

index  of  refraction  (these  fluctuations  being  caused  by  atmospheric  turbulence 

mixing  effects)  and  I  is  the  approximate  size  of  the  largest  eddies.  It  will  be 
2° 

explained  that  the  CR  parameter  can  only  be  employed  in  a  certain  domain  of 
turbulence  scales  (eddy  sizes),  namely  the  "inertial"  range  in  which  there  is  a 
"cascade"  of  energy  flow  from  the  large  eddy  scale  sizes  to  the  small  ones  where 
the  viscous  dissipation  occurs.  The  so-called  "structure  function"  of  the  index  of 

refraction,  r),  (defined  as  the  mean  square  difference  of  index  of  refraction 

n  2  2/^ 
measured  at  two  places  a  distance  r  apart)  is  equal  to  r  '  in  the  inertial  range. 

Section  2  displays  both  an  heuristic  and  a  mathematical  treatment  of  the  turbulence 

2 

concepts  which  make  up  the  context  of  C  (the  mathematical  treatment  being  given 
in  an  appendix). 

2 

Section  3  is  devoted  to  the  display  of  examples  of  how  Cn  affects  ootical  sys¬ 
tems.  The  turbulent  eddies  can  be  described  as  a  random  array  of  lenses  or  blobs 


3 


of  air  with  index  of  refraction  differing  slightly  from  their  surroundings.  It  is 
2 

shown  that  the  Cn  effects  on  electromagnetic  propagation  can  be  fairly  accurately 
explained  on  this  simple  basis.  The  two  main  effects  are  phase  modulations  (  and 
their  associated  wavefront  distortions),  and  "scintillation"  which  is  the  random 
fading  in  and  out  of  the  signal.  The  "twinkling"  of  stars  is  precisely  this  effect 
(also  known  to  astronomers  as  "bad  seeing"  conditions).  Scintillation  can  be  ex¬ 
plained  as  being  due  to  the  focusing  and  defocusing  effects  of  the  "blob  lenses" 
already  mentioned. 

Section  4  describes  the  current  state  of  the  art  of  estimating  C  by  means  of 

standard  "weather  balloon"  (that  is,  raw  insonde-radiosonde)  observations.  The 

objective  of  such  research  is  to  convert  the  entire  data-bank  of  such  observations 

2 

into  a  comprehensive  set  of  estimates  of  C  over  a  wide  range  of  spatial  locations 
and  times.  Thi3  would  make  possible,  presumably,  climatological  and  forecasting 
capabilities  for  C^. 

Section  5  considers  the  problems  to  be  solved  by  future  research,  and  possible 

techniques  to  consider.  In  essence,  the  weakest  part  of  the  current  analysis  for 

radiosonde  measurements  concerns  the  lowest  altitudes  (G  km  to  ground).  More 

experimental  investigation  is  urgently  needed  here.  Another  question,  the  most 

2 

crucial  in  fact,  is  whether  or  not  there  exists  a  significant  difference  in  C  as  a 
function  of  location  and  season.  "Anecdotally.  "  the  placement  of  a  telescope  site 
for  optimum  "seeing"  conditions  would  seem  to  indicate  that  C~  does  depend  on  lo¬ 
cation;  but,  at  this  writing  the  question  does  not  seem  to  be  answered  to  everyone's 

satisfaction.  In  contrast,  it  has  already  been  established  that  there  is  a  significant 

2  2 
diurnal  effect  on  C^.  For  the  purpose  of  making  low  altitude  measurements  by 

radar,  it  will  be  necessary  to  make  use  of  a  powerful  radar  system  equipped  with 
a  steerable  dish  (Lincoln  Labs,  for  example,  has  such  a  piece  of  equipment).  In 
addition,  the  "rnicroshears"  which  cause  the  C  -producing  turbulence  can  be  meas¬ 
ured  by  in-house  generated  technology  involving  smoke  trails.  The  latter  type  of 

2 

measurements  would  be  useful  for  verifying  current  C  models  for  all  altitudes. 
Balloon-borne  instrumentation  techniques  (also  generated  in-house)  should  also  be 
very  useful. 

The  report  concludes  with  an  annotated  bibliography  for  those  who  may  wish 
to  read  further  on  this  subject. 
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Optical  Turbulence  Forecasting: 
A  Tutorial 


1.  INTRODUCTION 

Astronomers  have  long  known  that  atmospheric  turbulence  can  cause  optical 
propagation  effects.  The  famous  "twinkling"  of  stars  is  due  to  such  turbulence. 
When  turbulence  is  at  a  minimum,  astronomers  can  photograph  planets  etc.,  with 
a  minimum  of  blur,  because  with  less  turbulence  there  is  less  so-called  "image 
dancing."  Brightness  fluctuations,  or  scintillations,  also  diminish  when  there  is 
low  turbulence.  In  general,  such  conditions  are  called  "good  seeing  conditions." 
The  basic  cause  of  the  optical  effects  is  the  presence  of  fluctuations  in  the  index 
of  refraction  along  the  optical  path.  These  fluctuations  are  caused  by  temperature 
fluctuations  and,  in  the  case  of  microwave  frequencies,  fluctuations  of  humidity. 
The  inhomogeneous  nature  of  such  quantities  is  due  to  the  (incomplete)  mixing 
effects  of  turbulence.  In  a  homogeneous  or  "already  totally  mixed"  region  of 
atmosphere,  turbulence  will  have  no  optical  effect.  So,  variations  of  index  of 
refraction  can  be  caused  by  turbulence;  but  there  can  be  (in  principle)  turbulence 
without  index  of  refraction  variations.  The  time  rate  of  optical  modulations  due  to 
the  fluctuations  depends  both  on  the  sizes  of  the  eddies  (or  inhomogeneities)  and  on 
the  mean  motion  of  air  which  carries  these  fluctuations  across  the  optical  path. 

The  Air  Force  is  currently  involved  with  communication  and  navigation  sys¬ 
tems  which  use  laser  beams  to  carry  signals.  Such  systems  can  involve 

(Received  for  publication  18  January  1980) 
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ground-to-space  and  space-to-ground  propagation  as  well  as  alr-to-space,  ground- 
to-ground,  etc.,  links  as  well.  The  capabilities  of  such  systems  could  be  signif¬ 
icantly  affected  by  the  turbulence -induced  fluctuations  mentioned  above.  In  addition, 
photography  from  high  altitudes  can  also  be  influenced  by  these  effects  as  well  as 
any  weapon  system  which  depends  upon  the  propagation  of  a  laser  beam.  Refrac¬ 
tion  fluctuations  are  not  bad  for  all  Air  Force  propagation  systems,  however.  For 
example,  "over  the  horizon"  or  "troposcatter"  communication  systems  cannot 
function  at  all  unless  the  fluctuations  are  present.  For  these  reasons  it  would  be 
desirable  to  have  some  way  to  forecast  those  conditions  which  give  rise  to  the 
fluctuations.  Two  important  and  interesting  questions  are  whether  there  are  geo¬ 
graphic,  diurnal,  seasonal,  or  weather  conditions  which  are  especially  conducive 
to  fluctuations  (or  to  their  absence)  and  whether  the  current,  large  data  base  of 
radiosonde  measurements  can  be  tied  to  the  value  of /C^  (a  measure  of  index  of 

refraction  fluctuations)  in  such  a  way  that  it  can  be  used'to  investigate  the  depend- 
2 

ence  of  Cn  upon  other  measurable  —  perhaps  predictable  —  parameters.  One  pur¬ 
pose  of  this  report  is  to  examine  the  current  state  of  the  art  along  these  lines.  In 

particular,  we  shall  examine  the  best  methods  available  to  relate  radiosonde  esti- 
2 

mates  of  C  as  a  function  of  altitude  and  time  to  simultaneous  colocated  radar 
n 

measurements. 

It  is  hoped  that  a  non-expert  in  this  field  will  be  able  to  see  from  this  report 
2 

how  C  is  related  to  turbulence-induced  index  of  refraction  fluctuations  and  how  the 
n 

latter  can  affect  system  performance.  Also  he  may  see,  hopefully,  how  one  may 
someday  succeed  in  connecting  rawinsonde  measurements  to  the  relative  perform¬ 
ance  of  such  systems. 


2.  C??  AND  TURBULENCE-INDUCED  INDEX  OF  REFRACTION 
FLUCTUATIONS 

In  Section  3  we  shall  consider  the  •'onnection  between  index  of  refraction  fluc- 

2 

tuations,  as  parameterized  by  C  ,  and  optical  system  performance.  It  will  be 
shown  that,  in  many  ways,  the  turbulence-induced  inhomogeneities  can  act  like  a 
random  assembly  of  lenses.  The  signal  amplitude  changes,  for  example,  can  be 
explained  in  terms  of  focusing  and  defocusing  effects  due  to  these  lenses  as  they 
are  convected  across  the  optical  path.  In  the  present  section,  however,  we  shall 
concentrate  on  the  physics  behind  the  parameter  C“  which  will  appear  in  all 
equations  describing  the  effects  of  turbulence  on  electromagnetic  propagation.  Our 
starting  point  is  the  concept  of  the  "structure  function," 
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2.1  The  Structure  Function 

The  index  of  refraction  depends  exclusively  upon  temperature  for  some  elec¬ 
tromagnetic  frequencies  (that  is,  for  optical  frequencies)  whereas  for  radar  fre¬ 
quencies  it  depends  on  both  temperature  and  specific  humidity.  For  simplicity  we 

shall  initially  ignore  this  double  dependence  and  Imagine  the  case  where  only  the 
* 

temperature  is  important  (the  dependence  upon  water  vapor  will  be  treated  at  the 
end  of  Section  2).  For  generality  we  can  start  by  considering  an  arbitrary  passive 
scalar  and  call  it  "0. "  The  structure  function  is  defined,  not  in  terms  of  0,  but  in 
terms  of  6',  that  is,  the  "fluctuation"  of  8,  which  is  defined  by 

8  a  7  +  O'  .  (1) 

Thus  0'  is  the  deviatior  from  the  mean  of  9  (denoted  by  9).  The  structure  function, 
for  9  is  thus 

£%(r)  a  >r)  -  fl'ir^]2)  .  (2) 

The  angular  brackets  {  )  signify  an  ensemble  average  but  in  practical  terms  it 
may  in  fact  be  a  time  average.  The  vector  T  is  the  distance  between  the  two 
locations  of  the  measurements:  "rj  +  7-  and  r^.  When  the  turbulence  is  isotropic 
and  homogeneous,  «0^(r)  will  depend  neither  upon  location  rj  or  the  direction  of 
F.  We  shall  consider  only  such  simplified  turbulence  from  now  on.  Fortunately 
this  is  a  good  approximation  for  our  purposes.  For  this  reason,  r  will  be  repre¬ 
sented  as  a  scalar. 

During  the  development  of  turbulence  theory  by  such  innovators  as  G.  I.  Taylor, 
work  was  done  with  the  so-called  autocorrelation  function,  B0(r),  which  is  defined 
by 


B0(r)  s  <fl'(r1  +  r)  e'(r,»  . 

The  quantities  B.(r),  and  <5^(r)  are  closely  related. 

U  V 

J  «^<r)  -  (B0(o)  -  B0(r)l  . 


(3) 

In  particular  (Tatarski,  *  p.  10) 

(4) 


1,  Tatarski,  V.  I.  (19GI)  Wave  Propagation  in  a  Turbulent  Medium.  McGraw- 
Hill,  N.Y.  * 

♦ 

Temperature  in  the  following  will  always  be  considered  as  a  passive-scalar  addi¬ 
tive  to  the  fluid  even  though  in  actual  fact  it  is  not.  The  approximation,  however, 
in  the  cases  of  present  interest,  is  a  good  one. 
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The  structure  function,  a  Russian  development,  has  now  acquired  universal  usage 
because  of  certain  practical  advantages  it  affords.  Whereas  BJr)  is  advantageously 
used  in  connection  with  turbulence  theory  (due  to  its  very  simple  connection  with 
the  spatial  frequence  power  spectra  which  are  intrinsic  to  much  of  that  theory)  the 
structure  function  has  the  following  experimental  advantages  <Panchev.2p.  69).  The 
fluctuation  measurements  are  almost  always  associated  with  "slow  variations"  or 
’trends."  The  structure  constant,  being  a  first  difference,  automatically  filters  out 
hese  annoyances. "  It  is  also  less  sensitive  to  error  due  to  random  deviations  It 
also  provides  advantages  for  instrumental  design  in  turbulence  measurements. 
Finally,  ^(r)  has  the  advantage  that,  theoretically,  it  has  a  beautifully  simple 
form  when  one  is  dealing  with  "inertial  range"  or  "cascade"  turbulence  It  is 
precisely  this  simplicity  which  gave  rise  to  the  concept  of  C2.  We  shall  in  the 
present  context  call  it  C2.  and  this  more  general  term  is  called  the  structure  con¬ 
stant  (or  rather  Cfl>  not  squared,  is  the  structure  constant).  As  will  be  shown 
below  for  inertial  range  turbulence, 

•%<r>  »  C2  r2/3  , 

e  e  • 


Equation  (5)  was  the  original  form  given  by  Kolmogorov;  however,  the  one-dimen- 
sional  spectrum  which  is  associated  with  Eq.  (5)  depends  on  the  spatial  frequency 
Mor  (2n) /wavelength)  as  k  One  usually  hears  of  "Kolmogorovs  minus  five- 

thirds  law  in  truth  he  discovered  the  "two-thirds  law"  of  Eq.  (5). 

The  relationship  between  <%(r>  and  Bp(r)  as  well  as  the  variance,  a2  is 
graphically  illustrated  in  Figure  1  (cf.  Panchev,  2  p.  69).  9 

Now  that  &£r)  has  been  defined  we  can  now  travel  on  a  path  through  turbu¬ 
lence  theory  which  leads  to  a  rather  physical  derivation  of  Eq.  (5)  (see  Appendix  A 
for  more  detail).  Only  "cascade"  turbulence  will  be  considered 


2.  P^chev,  S.  (1971)  Random_FunctionS  and  Turbulence  Pergamon  Press, 


2.2  Tha  Inertial  Range  amt  JEhr)  oc  C*  t2** 

In  the  following  we  shall  consider  a  simplified  form  of  turbulence  theory.  The 

main  objectives  are  to  derive  the  structure  function  from  cascade  theory  and  to 

relate  the  structure  constant  to  the  energy  dissipation  rate  of  the  turbulence  ,  e, 

and  the  "dissipation  rate,"  related  to  the  molecular  smoothing  out  of  Inhomo- 

M  2 

geneltles  so  that,  at  the  end,  v»«s  can  arrive  at  the  relation  for  Cn  used  by 

VanZandt  et  al^  In  their  work  on  radiosonde  estimation  of  cj!.  Our  treatment 

1  0 

mathematically  follows  Tatnrakt  but  the  physical  discussions  are  added. 

For  the  sake  of  definiteness  we  shall  replace  the  meaning  of  8,  which  has 
been  "passive  scaler"  up  to  this  point,  with  potential  temperature."^  In  this  way  9 
will,  in  the  case  where  humidity  la  not  Important,  be  closely  related  to  Index  of 
refraction.  The  exact  connection  ia  given  at  the  end  of  Section  2, 

Next,  we  make  certain  simplifying  assumptions.  In  addition  to  homogeneity 
and  Isotropy  of  turbulence  mentioned  above,  we  Introduce  "steady  state,"  that  Is, 
we  assume  that  as  much  energy  flows  Into  tne  turbulence  as  13  dissipated  In  the 
form  of  heat  (due  to  viscous  dissipation).  Furthermore  we  shall  resort  to  the 
famous  "cascade"  picture,  first  introduced  by  Richardson  in  the  form  of  a  poem 
and  then  put  effectively  into  mathematical  form  by  Kolmogorov  (and  subsequently 
but  independently  by  Onsager  and  von  Weizsflcker).  In  the  so-called  "inertial 
range"  the  large  eddies  break  down  into  smaller  eddies,  due  to  fluid  dynamic 
instabilities,  and  these  smaller  eddies  likewise  break  down  Into  still  cmaller  ones. 
This  "cascade  of  energy"  from  large  sizes  or  "scales"  to  smaller  scales  continues 
until  the  scale  becomes  so  small  that  viscosity  effects  become  important,  at  which 
point  the  energy  flows  out  of  the  "cascade"  and  Into  the  ''sink"  of  molecular  heat. 
As  Richardson  said: 

"Big  whirls  have  little  whirls. 

Which  feed  on  their  velocity; 

And  little  whirls  have  lesBer  whirls, 

And  so  on  to  viscosity  (in  the  molecular  sense)," 

The  simplest  derivation  of  the  structure  function  in  this  inertial  range  can 
now  be  given  provided  one  knows  that,  in  addition  to  r,  only  two  other  parameters 
are  of  any  physical  relevance,  namely  e  and  eN  mentioned  above.  This  is  true 


e  represents  heating  due  to  turbulence  and  represents  molecular  smoothing 
(diffusion)  of  turbulent  inhomogeneities. 

^Potential  temperature  is  used  instead  of  temperature  in  a  compressible  gas  in 
order  to  remove  the  complication  that  a  drop  in  pressure  decreases  the  tempera¬ 
ture;  the  potential  temperature  in  this  case  would  remain  constant. 

3.  VanZandt,  T.  E. ,  Green,  J.L.,  Gage,  K.  S. ,  and  Clark,  W.  L.  (1978)  Vertical 
profiles  of  refractlvity  turbulence  structure  constant:  Comparison  of 
observations  by  the  Sunset  Radar  with  a  new  theoretical  model.  Radio 
Science  13:818-829. 


because  we  have  a  steady  state  and  the  cascade  must  depend  only  on  the  rates  of 
dissipation  (heat  =  e,  smoothing  =  eN)  at  the  small  end.  Once  this  Is  granted, 
one  can  write 


&(r)  .  f!N  r2/3 
S'  '  el/3 


(6) 


This  follows  from  dimensional  analysis.  The  symbol  a  is  for  a  constant  of  order 
unity.  We  shall  start  again,  however,  from  the  "equations  of  motion"  and  rederlve 
Eq.  (6)  with  a  little  more  detail  In  order  to  gain  more  insight  into  its  meaning  and 
to  obtain  information  which  will  enable  us  to  put  Eq,  (6)  into  a  form  more  ammen- 
able  to  experimental  requirements. 

The  starting  point  is  the  molecular  diffusion  equation.  This  is,  in  fact,  the 
"equation  of  motion"  for  the  passive  scalar  (potential  temperature  in  our  case),  0. 

gf  +  dlv  (-Vq  grad  9)  =  0  ,  (7) 


where  is  the  molecular  coefficient  of  diffusion  for  6,  and  where 


d9 

3f 


+  v  •  grad  9 


(8) 


that  is,  the  "mobile  derivative"  and  v  the  velocity,  div  is  the  divergence,  and  grad 
the  gradient.  Since  it  turns  out  that  the  dynamics  will  not  be  significantly  affected 
by  compressibility,  we  make  the  usual  assumption  that  the  air  can  be  considered 
incompressible  in  this  context.  Mathematically,  then,  we  put 

div  (v)  =  0  (9) 

and  this  makes  it  possible  to  write 

v  •  grad  (0)  =  div  (v0)  .  (10) 

Thus  Eq.  (7)  becomes 

|f  +  div  (v0  -  v6  grad  9)  =  0  (11) 


and  this  is  the  basic  equation.  At  this  point  we  again  consider  fluctuations. 
Reynolds  was  the  first  to  introduce  the  "decomposition"  of  a  variable  into  its  mean 
and  fluctuating  part.  We  do  this  for  both  9  and  v^  where  v.  is  the  i^  component  of 
the  velocity  vector,  v. 
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jra^atraE^^r"""'  ^■TTB^,vy.,aaiwnria«»waai 


6  m  6+0' 

(12) 

vt  .  vt  +  v| 

whore,  as  before,  the  overbar  is  an  average.  In  Appendix  A  It  is  shown  how  the 
insertion  of  Eq.  (12)  Into  Eq.  (11)  and  subsequent  multiplication  by  O'  and  aver- 
aging  etc.  lead  to  a  relation  between  (O'  )  (that  is,  the  mean  square  fluctuations 
of  6')  and  e^.  The  relation 

Ke  (grad  e)2  =  vQ  ((grad  0’)2>  (13) 

is  the  result  from  such  manipulations.  The  angular  brackets  are  an  average,  and 

Kg  is  the  so-called  eddy  diffusivity  to  be  discussed  and  defined  below.  The  right 

side  of  Eq.  (13)  represents  the  rate  that  the  average  fluctuation  or  inhomogeneity, 
— 

6'  ,  is  wiped  oxst  per  unit  of  time  by  molecular  diffusivity,  v  g.  The  right  side  de¬ 
fines  e^.  The  left  side,  which  we  shall  call  P,  represents  the  "production  term," 
that  is,  the  rate  of  production  of  the  fluctuations  per  unit  of  time. 

P  =  Kg  (grad  0)2  .  (14) 

The  eddy  diffusivity,  K_,  can  be  approached  in  several  ways.  Its  dimensions 

9-1  ^ 

are  (L  ][T  M.  It  could  be  written  as 


Kg  =  v'l 1  (15) 

that  is,  as  a  product  of  a  turbulent  velocity  scale  and  a  turbulent  length  scale. 

4 

This  is  its  basic  physical  meaning.  However,  this  leaves  the  value  of  t '  unknown 
and  for  present  purposes  we  can  solve  this  problem  by  defining  Kg  by 

Kg  -  -  (0'  v1)  /(grad  6)  .  (10) 

While  it  is  clear  that  the  dimensions  of  Kg  are  correct,  (v1  •  /'),  and  that  t '  is 
appropriately  defined  in  terms  of  a  length  relating  to  0  variation,  it  is  not  obvious 
yet  how  one  arrives  at  Eq.  (10)  directly.  The  best  way  to  derive  Eq.  (10)  is  to 
make  use  of  an  analogy  with  Fourier's  law  Tor  heat  flow.  Calling  the  flux  of  tem¬ 
perature  flow  q  we  have 


4.  Pasquill,  F.  (1974)  Atmospheric  Diffusion  (2nd  Ed.),  John  Wiley  &  Sons, 
New  York. 
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qfl  -Kfl  grac!  0  ■  (17) 

Rut  the  turbulent  flux  qff  la  just  the  flow  of  the  fluctuations  of  O'  aa  carried  by  the 
velocity  fluctuations  v1.  This  in  turn  is  given  by 

q0  •■=  <v'  O')  (18) 

which  physically  means  that  the  flux  is  due  to  the  correlation  between  velocity  and 
potential  temperature  fluctuations.  Substitution  of  Eq.  (18)  into  Eq.  (17)  (that  is. 
eliminating  q0)  one  arrives  at  Eq.  (16)  directly. 

At  this  point  we  wish  to  note  a  certain  analogy.  Eq.  (14)  defines  the  produc¬ 
tivity  of  turbulent  fluctuations  while,  as  we  have  already  stated, 

eN  ■  vd  «grad  0')2>  .  (19) 

Notice  the  similarity  in  form  for  P  and  eN-  We  shall  make  use  of  this  shortly. 

Next  we  shall  consider  the  physical  meaning  of  P.  After  all,  we  did  not  derive 
Eq.  (13)  here  explicitly  (see  Appendix  A),  therefore  we  resort  to  "understanding" 
it  from  its  physical  interpretation.  To  do  this  we  start  by  making  a  simplification 
which  turns  out  to  be  the  usual  and  justifiable  one  for  the  atmosphere,  namely  that 
the  mean  value  of  6,  that  is,  0,  depends  only  upon  the  altitude,  z.  In  this  way  we 
can  write 

p ■  <e'  v;,>  f  ■  (so 


Now  we  shall  examine  why  P  is  physically  the  rate  of  production  of  6'2  per  unit 
time.  Figure  2  makes  possible  a  very  simple  explanation  of  this.  In  this  figure, 
we  imagine  a  parcel  of  air  originally  located  at  altitude  Zj  and  let  it  be  transported 
(by  an  eddy,  say)  to  a  higher  altitude  z^  which  is  at  a  height  Az  above  z  .  We  shall 
assume  that  this  parcel  originally  had  the  average  temperature  b  at  z^.  At  its 
destination  z it  still  has  this  potential  temperature  (because  potential  tempera¬ 
ture  does  not  change  as  the  parcel  expands  to  a  larger  volume  at  the  higher  altitude 
and  because  we  assume  that  the  parcel  keeps  its  identity,  that  is,  it  does  not  mix). 
When  the  parcel  is  at  z2,  however,  its  temperature  now  differs  from  the  mean.  6, 
at  thi3  new  altitude.  This  difference  is  given  in  Figure  2  and  we  call  it  AO.  AO  is, 
in  effect,  the  value  of  O'  for  the  parcel.  Notice  the  important  fact  that  the  grad¬ 
ient,  d0/dz,  is  necessary  for  the  development  of  AO .  Mathematically 
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Figure  2.  Generation  of  Fluctuation,  O',  By 
Vertical  Transport 


(21) 


Now  multiply  Eq.  (21)  by  O'  and  perform  an  average,  (  ),  to  obtain 

{O'2)  =  (O'  Az)  £  .  (22) 

o 

In  order  to  estimate  the  rate  of  change  of  (O'  )  with  respect  to  time,  we  divide 
Eq.  (22)  by  At  and  identify  Az/At  with  the  z  component  of  the  eddy  velocity  fluc¬ 
tuations,  v' ,  to  obtain 
z 


%r  -  <0’  vk>  § 


and  this  is  precisely  Eq.  (20).  Using  the  definition  for  K0  given  by  Eq.  (16),  we 
arrive  at  Eq.  (14). 

It  is  worthwhile  to  restate  the  physical  meaning  of  P.  It  represents  a  con¬ 
version  of  the  gradient  of  0  (a  large  scale  variation)  into  the  fluctuations  O'  (a 
small  scale  variation)  of  the  turbulence.  It  is  brought  about  in  part  by  the  turbulent 
fluctuations  in  velocity  v'.  At  the  smallest  end  of  the  cascade,  the  fluctuations 
(O'  )  are  smoothed  out  by  molecular  diffusion.  The  latter  occurs  at  the  rate  e^. 
The  assumption  of  "steady  state"  is  nothing  other  than  the  requirement  given  by 

P  =  eN  .  (24) 
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Now  consider  Eq.  (19)  and  let  us  compare  It  to  P.  Compare,  in  other  words, 


(25) 


w  1th 


eN  =  ve  ^grad  e')2)  ■  <2r>) 

These  are  of  exactly  the  same  form.  In  a  manner  similar  to  P,  eN  can  be  viewed 
as  a  transition  between  a  relatively  large  scale  of  variation  of  8  to  a  small  scale. 
The  big  difference  is  that  relates  to  the  conversion  to  molecular  variation. 

Thus  K.  parameterizes  the  conversion  from  mean  variation  to  eddy  variation  and 
Vq  does  the  same  for  conversion  from  eddy  fluctuations  to  molecular  fluctuations. 
At  this  point  the  physical  meaning  of  Eq.  (13)  should  be  in  evidence. 

We  now  turn  to  the  next  step  in  deriving  the  structure  function.  This  consists 
on  making  an  estimate  of  (O'2)  in  terms  of  e^,  e,  and  the  scale  of  turbulence  or 
size  of  the  eddies  which  we  shall  here  call  i.  To  start  this  we  make  a  simple 
estimate  of  the  quantity  To  do  this,  consider  Eq.  (2G)  and  approximate 

<(grad  9')2)  by  {<e|2)/(i)2).  Thus 


e  N  "  f2 


(27) 


Physically  we  imagine  O'  as  the  fluctuation  of  0  across  the  eddy  of  size  /.  Quan¬ 
tity  t  here  is  defined  as  the  "representative  scale"  of  the  turbulence.  In  order  to 
have  a  cascade  of  the  type  discussed  above,  however,  the  molecular  dissipation 
rate  of  the  largest  eddies  must  be  very  much  smaller  than  the  rate  that  they  give 
up  their  fluctuations  to  the  next  smaller  size  of  eddies  ("big  blobs  become  smaller 
blobs").  Consider  an  eddy  of  scale  I.  It  converts  itself  into  a  group  of  smaller  blobs 
in  a  time  t^ .  It  is  known  in  turbulence  theory  that  this  time  t^  is  very  short.  In 
fact,  it  is  one  "turn  over  time"  of  the  eddy.  Let  be  the  internal  velocity  assoc¬ 
iated  with  the  eddy  (v^  is  the  same  thing  as  the  velocity  fluctuation,  v',  across  I). 
Then 


I 


Of  course,  e  must  be  interpreted  as  a  "passive  scalar"  in  this  context  since  tem¬ 
perature  is  meaningless  at  the  molecular  level. 


Now  in  order  that  a  cascade  exist.  It  Is  required  that 


«  l?i>  - 


‘N - 1 


(29) 


but  from  Eq.  (27)  and  Eq.  (15)  we  see  that  Eq.  (29),  written  as 


'N 


■  v  «  V  i  m  ^ 

I 2  6  |2  *  1 2 


(30) 


which  Implies  that  for  a  cascade  the  condition  Is 


Ke  >>  ve 


(31) 


In  other  words,  the  big  eddies  must  break  up  at  a  much  faster  rate  than  the  action 
of  molecular  effects  upon  them.  Alternatively,  diffusion  due  to  turbulence  is  much 
greater  than  that  due  to  molecular  transport.  At  the  small  scale  end  of  the  cas¬ 
cade,  where  the  scales  are  so  small  that  molecular  effects  are  relatively  large, 
we  have  the  condition 


where  v^  and  1  j  are  the  velocity  and  spatial  scales  of  the  dissipating  small  scale 
eddies.  This  defines  the  "dissipation  scale,"  ij  of  the  cascade. 

We  need  another  relationship  for  eN  because  we  wish  to  relate  it  to  e.  To  do 
this  we  use  the  steady  state  condition,  Eq.  (24),  or 


*N  ■  P  '  <»’  '”>  § 

and  make  the  approximation  that  (9'  ~  l  (d?/dz)).  In  this  way 

{O'2)  v 
6  N  t 


(33) 


(34) 


Now,  recall  the  definition  of  e.  It  is  the  rate  that  the  kinetic  energy  per  unit  mass 
of  the  velocity  fluctuations  (v1  )  becomes  converted  into  heat  through  molecular 
dissipation  effects.  We  need  an  estimate  for  e  in  terms  of  v  ^  in  order  to  eliminate 
the  latter  from  Eq.  (34).  The  physical-dimensional  argument  goes  as  follows. 

The  time  scale  for  the  conversion  of  an  eddy  into  a  number  of  smaller  ones,  that 
is,  the  time  scale  for  an  eddy  to  give  up  its  energy,  is  (t /v()  as  has  already  been 
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explained,  Suppose  now  that  all  of  the  kinetic  energy  of  the  large  eddies  of  sice  I 
Is  transfered  down  to  the  next  smaller  size  of  eddies  In  this  time;  then,  this  would 
represent  (in  steady  state)  the  rate  of  energy  flow  (per  unit  mass)  into  the  top  end 
of  the  cascade  of  (v^)  per  Interval  of  time  (i /v^).  Ultimately,  In  analogy  with 
P  *  eN,  this  must  equal  e.  Thus, 


e 


vf 

T 


(35) 


This  equation  for  e,  In  terms  of  the  scales  v  and  I,  as  simple  as  It  may 

I  » 

appear  and  as  "dimensional"  and  non-rigorous  as  It  may  seem  at  first  sight,  Is 

5 

one  of  the  "corner  stone"  assumptions  of  turbulence  theory.  It  la  due  to  the 
validity  of  Eq.  (35)  that  one  can  say  turbulence  is  a  "strongly  damped,  non-linear 
stochastic  system."  The  dimensional  arguments  which  are  characteristic  of  our 
treatment  up  to  this  point  are  also  typical  of  much  of  turbulence  theory  (due  to  the 
difficulties  Introduced  by  non-linearity).  If  they  appear  repugnant  to  the  reader 
who  may  prefer  exact  analytic  approaches,  it  should  be  pointed  out  that  it  is  pre¬ 
cisely  these  "similarity"  techniques  which  can  claim  a  large  part  of  the  useful 
practical  aspects  of  turbulence  theory  (taken,  of  course,  in  conjunction  with 
experimental  measurements),  t 

Now  using  Eq.  (35)  to  eliminate  from  Eq.  (34) 

<e'2)  (ei)1/3 

=  j  (36) 

2 

and  solving  for  (ej  )  we  have 


e  f2/3 
eN* 


J7 T 


(37) 


Eq.  (37)  is,  in  effect,  an  estimation  of  the  structure  function  for  the  special 
case  where  r  -  t.  For  smaller  r,  the  structure  function  must  be 

★ 

Vj  and  i  j  refer  to  dissipation  scales  to  be  explained. 

^The  analytic  approach  has  been  worked  out  for  the  inertial  range.  See  .Appendix  B. 

5.  Tennekes,  H. ,  and  Lumley,  J.L.  (1972)  A  First  Course  in  Turbulence, 

MIT  Press,  Cambridge,  MA. 


18 


(38) 


f 

i 


( 

I 

!) 


r2/3 


where  a  Is  a  constant  of  order  unity.  The  range  of  validity  of  Eq.  (38)  is  confined 
to  the  cascade  range.  Denoting  l  as  the  size  of  the  largest  eddies  or  "outer 
length,  "  or  the  size  of  the  "energy-containing  eddies,  "  or  the  so-called  "integral 
scale",  and  denoting  the  small-scale  end  of  the  cascade  by  l  ^  (cf.  Eq.  (32)),  we  can 
denote  the  range  of  validity  of  Eq.  (38)  by 


(1.  «r«l  ) 
l  o 

and  using  v(  from  Eq.  (35)  in  Eq.  (32)  we  can  solve  for  I  j  to  obtain 


This  length,  <j,  is  called  the  "inner  length"  or  "microscale." 
Eq.  (38)  is  also  written  as 


o£^(r)  =  C2  r2/3  (40) 

where  Cfl  is  called  the  structure  constant  as  mentioned  in  Section  2,  thus  we  have 
the  relation  (see  Eq.  (38)) 


ae. 


.1/3 


(41) 


2  2 
The  quantity  Cg  is  an  object  which  very  closely  resembles  the  quantity  CR  as  we 

shall  see. 

2 

What  is  the  physical  interpretation  of  C_?  From  Eqs.  (37)  and  (41)  we  find 

U 


2  (0'2> 
0  72T 


(42) 


fi.  Hinze,  J.O.  (1975)  Turbulence  (2nd  Ed.),  McGraw  Hill. 
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Thus,  Cg  la  the  ratio  of  the  mean  square  fluctuations  of  6  to  the  2/3  power  of  the 
dimension  t  of  the  eddies.  In  Eq.  (42)  the  term  I  means  the  representative  length 


scale  of  the  eddies.  Thus  we  can  also  write 


«2£<r)  -  <e<2>  (43) 

f) 

where  ( O' **)  corresponds  to  the  variations  over  a  distance  1  of  the  scale  of  turbu¬ 
lence.  When  r  =  I,  =  (0|2). 

For  completeness  we  add  the  form  for  ,£^(r)  which  is  valid  for  the  case  where 
r  «  f  j.  It  is  (Tatarskl,  *  p.  46) 

■  <«> 

In  order  to  Illustrate  the  experimental  significance  of  the  structure  function 
and  the  structure  constant  we  consider  a  specific  example.  Consider  a  device 
which  measures  temperature  at  two  points  separated  by  a  fixed  horizontal  distance 
r1.  (We  shall  let  0  represent  ordinary  temperature.)  Assume  that  r'  is  chosen  so 
that  it  satisfies  /  ^  «  r'  «  l  for  the  case  of  interest.  Suppose  further  that  this 
device  (a  pair  of  temperature-measuring  instruments  in  effect)  is  attached  to  a 
balloon  and  sent  aloft  to  measure  the  ,j^(r)  of  various  altitude  regions  of  the 
atmosphere.  The  first  question  we  wish  to  answer  is:  "How  is  the  average  in 
Eq.  (2)  actually  performed,  the  latter  equation  giving  the  definition  ofo&'lr)?" 

More  specifically  we  wish  to  know  whether  it  is  a  spatial  average  or  an  average 
over  time.  The  answer  is  that  it  is  a  time  average;  but,  since  in  general  such 
measurements  must  be  made  from  a  location  relatively  far  below  the  balloon  (to 
avoid  balloon-induced  artifacts),  there  will  be  a  relative  wind  velocity  between  the 
measuring  device  and  the  air  which  immediately  surrounds  it.  This  would  be  due 
to  the  fact  that  in  general  there  are  wind  shears  present  and  while  the  balloon 
might  be  nearly  at  rest  with  respect  to  its  surroundings,  the  gondola  below  would 
be  expected  to  be  moving  with  respect  to  its  environment.  Very  often  in  such 
situations  the  so-called  "frozen  turbulence"  hypothesis  is  used.  This  means  that 
the  measuring  device  traverses  the  turbulent  field  with  enough  velocity  to  render 
the  turbulent  velocity  fluctuations  very  small  by  comparison.  Thus  one  can  imagine 
that  the  turbulent  eddies  remain  fixed  or  "frozen  in  place"  as  the  measuring  device 
traverses  them.  Under  these  conditions  the  time  average  in  the  definition  of 
«£^(r)  is  actually  a  spatial  average  so  far  as  its  physical  significance  is  concerned. 

The  quantity  (S'2)  is,  ideally,  a  one-point  ensemble  average.  It  is  given  by 
B„(o).  From  Eq.  (4),  however,  we  see  that  this  implies  that 


•^(r)  =  <0'2)  ~  Be(r) 


(45) 


In  view  of  this,  why  physically  ia  j£^(f)  =  (0'2)?  The  answer  is  that  when  r  reaches 

beyond  the  value  of  the  representative  scale  f,  the  correlation  Bg(r)  falls  off 

rapidly  (this  in  fact  is  related  to  the  empirical  definition  off).  Thus,  to  a  good 

approximation  we  can  imagine  that  when  r  is  somewhat  larger  than  f,  Eq.  (43)  no 

longer  applies  (we're  out  of  the  inertial  range)  and  B^(r)  goes  to  zero  because  there 

is  no  correlation  when  r  is  larger  than  the  existing  energy-containing  eddies.  This 

2 

completes  our  discussion  of  the  physical  significance  of  C  .  We  now  turn  to  the 

2  n 

problem  of  converting  C.  into  a  form  more  amenable  to  experimental  measure- 

(7 

ment. 

2.3  More  Useful  Forms  for  Calculations 

From  Eqs.  (24)  and  (25)  we  have 


One  can  derive,  in  a  very  similar  manner,  the  following  relation  for  e 


E  "  Si 


(47) 


where  U  is  the  mean  horizontal  velocity  and  K-.  is  the  eddy  diffusivity  for  momen- 
5  7  8  9 

turn  transfer.  " •  ’  •  By  inserting  these  forms  for  e  and  into  the  equation  for 
2 

Cg  Eq.  (41),  we  obtain 


7,  Lumley,  J.L.  ,  and  Panofsky,  H.A.  (1964)  The  Structure  of  Atmospheric 

Turbulence,  Interscience  Pub.,  N.  Y. 

8,  Lin,  J.T.,  Panchev,  S. ,  and  Cermak,  J.  (1969)  Turbulence  spectra  in  the 

buoyancy  subrange  of  thermally  stratified  shear  flows.  Project  Themis. 
Tech.  Report  No.  1,  College  of  Engineering,  Colorado  State  llniv. 

9,  Dewan,  E.M.  (1976)  Theoretical  Explanation  of  Spectral  Slopes  in  Strato- 

spheric  Turbulence- bata  and  Implications  for  Vertical  Transport. 
AFGL-TR -76-0247,  ADA036  307. 
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As  we  have  already  seen,  eddy  diffusivity  is  essentially  the  product  of  the  velocity 
and  spatial  scale,  that  is,  v^  /.  Thus  we  can  use 


Si  '  <vf 


(49) 


A 

(where  subscript  "o"  refers  to  outer  scale)  for  the  momentum  diffusion  coefficient. 
The  velocity  scale  for  momentum  transfer  is  given  by 


Eq.  (50)  is,  in  effect,  the  definition  of  the  velocity  scale  (or  the  "velocity  of  the 
energy  containing  eddies").  Thus 


Solving  this  for  (dU/dz)  and  inserting  into  Eq.  (48) 


(51) 


(52) 


2  2 

This  equation  for  is  the  form  used  in  estimating  by  means  of  radiosondes 
to  be  described  below. 

9 

We  now  turn  to  the  actual  value  of  C  .  The  latter  is  given  by  Eq.  (52)  when  8 
is  the  index  of  refraction.  The  latter  (in  the  case  of  radar)  depends  on  both  poten¬ 
tial  temperature  and  specific  humidity.  The  latter  dependence  is  especially  strong 
at  the  lowest  altitudes.  In  Eq.  (52)  (de/dz)  must  be  replaced  by 


where  N  is  defined  as  (n  -  1)  10^  and  where  n  is  the  index  of  refraction.  Here  0 
is  the  average  potential  temperature  and  q  is  the  ratio  of  the  mass  of  water  vapor 
per  unit  volume  to  the  mass  of  the  moist  air  per  unit  volume.  From  such  consid¬ 
erations,  Tatarski*  (p.  57)  has  shown  that 


In  shear  layer  turbulence,  f  *  i  but  /  >  1 .  There  seems  to  be  some  vagueness 

in  the  definition  of  "outer  length”"  Hef*e  l  would  be  layer  thickness  and 

'  3  (vrms)/(dU/dz)- 
-  2.8  (Ref.  3). 
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a 


(54) 


c;  =  a  (  P-  )ti  4/3  M2] 
n  1%  ° 


where 


M»  ~79  x  l0'G  P  (i  l.SSXtO^  /d T 

t2  l  'r 


3T  +  a 


7800  dq/dz 


.  ,  /  1.  55  X  104  q' 

1  +  (  - T - 


(55) 


(see  Ref.  1,  p.  57)  where  is  the  adiabatic  lapse  rate  of  9.  8  K  per  1000  m,  T  is 
the  absolute  temperature,  and  P  is  the  atmospheric  pressure  given  in  millibars. 

o 

VanZandt  et  al  have  modified  this  slightly  to  make  it  more  convenient  and  have 
obtained 


M  =  -77.  6  X  10‘(>  (  ) 


-6  /  P  \  la  in  6 
9z 


9  t  n  q 

1  +  15500  q  |  !  .i  V  9* 


9  in  6 

9z 


(S6) 


where  6  is  the  average  potential  temperature.  The  natural  logarithm  is  designated 
by  t  n. 


3.  C.2  AND  ATMOSPHERIC  TURBULENCE  EFFECTS  UPON 
OPTICAL  SYSTEMS 


3.1  A  List  of  Optical  Effect* 

Since  the  laser  systems  under  consideration  all  involve  propagation  through  a 
portion  of  the  atmosphere,  these  would  all  be  subject  to  a  number  of  degrading 
influences.  In  Pratt ^  many  designs  of  receivers  are  described.  In  some  of  these 
receivers,  such  as  the  heterodyne  and  homodyne  types,  a  "photomixing  process" 
is  used  to  mix  the  received  laser  beam  with  a  locally  generated  laser  beam  in 
order  to  create  a  low -beat  frequency  in  the  process  of  signal  detection.  Such 
systems  would  be  very  sensitive  to  any  corruption  of  the  signal  phase.  On  the 
other  hand,  all  receivers  would  be  subjected  to  the  "fading"  or  scintillation  effects. 
The  following  is  a  somewhat  comprehensive  list  of  turbulent  optical  (or  electro¬ 
magnetic  in  general)  effects. 

1.  Beam  steering:  The  laser  beam  can  be  deviated  from  the  line  of  sight  so 
that  part  of  it  (or,  in  the  worst  case,  all  of  it)  will  miss  the  receiving  aperture. 


10.  Pratt.  W.  K.  ( 1969)  Laser  Communication  Systems,  John  Wiley  &  Sons,  Inc., 
N.  Y. 


This  is  most  likely  for  the  case  where  the  receiver  is  in  space  and  the  transmitter 
is  on  the  ground  (partly  because  a  deviated  beam  will  drift  farther  from  the  target 
in  proportion  to  its  subsequent  distance  of  travel).  Beam  steering  effects  arise, 
when  the  size,  i,  of  the  atmospheric  inhomogeneity  (what  we  have  called  the  eddy) 
is  larger  than  the  width  of  the  beam,  This  occurs  almost  exclusively  on  the  uplink 
of  an  earth  —  space  channel. 

2.  Image  dancing:  The  atmosphere  can  cause  a  modulation  of  the  angle  of 
arrival  of  the  beam's  wavefront.  This  will  cause  the  image  of  the  beam  in  the 
receiver  (or  on  the  photographic  plate  in  the  case  of  photography)  to  be  focused  at 
different  places  in  the  focal  plane.  In  a  photographic  image  this  would  cause 
blurring,  for  example.  In  a  communication  receiver  this  effect  would  necessitate 
a  larger  aperture  and  hence  introduce  more  noise. 

3.  Beam  spreading:  In  this  case,  numerous  inhomogeneities  distributed 
across  the  beam  cross-section  cause  many  small  angle  scattering  events  to  occur. 
This  has  the  effect  of  spreading  the  beam  energy  over  a  wider  cross-section. 

This,  in  turn,  dilutes  the  signal  at  the  receiver  and,  in  consequence,  introduces  a 
decrease  in  the  signal-to-noise  ratio.  In  this  case  the  beam  width  must  be  signif¬ 
icantly  larger  than  the  eddies,  therefore  beam  spreading  is  a  problem  for  the 
downlink  of  an  earth-space  laser  system  (receiver  on  the  ground  and  transmitter 
in  space  or  high  altitude  aircraft). 

4.  Spatial  coherence  degradation:  Losses  of  phase  coherence,  across  the 
wavefront  of  the  beam  occurs  due  to  the  inhomogeneities  in  the  beam's  path.  This 
rapid  change  of  phase  with  respect  to  the  radial  position  within  the  beam  cross- 
section  interferes  with  the  photomixing  process  mentioned  above  in  connection 
with  the  mixing  of  received  signal  and  local  laser  beams  for  detection  purposes. 

In  the  case  of  photographic  systems,  this  effort  can  cause  blurring. 

5.  Scintillation:  Within  the  beam  cross-section,  interference  effects  can 
cause  destructive  and  constructive  reinforcement.  This  causes  the  power  to  vary 
widely  from  point  to  point  within  the  cross-section,  and  it  is  a  spatial  form  of 
scintillation.  A  temporal  form  of  scintillation  also  exists  which  consists  of  a 
fading  in  and  out  of  the  signal  in  time.  These  effects  are  due  to  the  focusing  and 
defocusing  effects  of  the  turbulent  eddies.  As  will  be  explained  below,  scintillation 
or  amplitude  variation  effects  can  be  explained  by  regarding  the  eddies,  or  index 
of  refraction  variations,  as  a  random  distribution  of  weak  lenses  distributed  within 
a  spatial  volume.  This  volume  is  to  be  envisioned  as  being  converted  across  the 
beam  and  thus  causing  effects  that  depend  on  time. 

ti.  Temporal  phase  fluctuations:  These  can  introduce  a  spurious  modulation 
which  can  destroy  the  amplitude  modulation  (AIU)  characteristics  of  a  laser  or 
radio  signal. 
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The  above  list  illustrates  the  corrupting  effects  of  turbulence  on  electromag¬ 
netic  propagation.  Of  course,  in  scatter  communication  systems,  the  lack  of 

turbulence  effects  would  have  the  result  of  destroying  the  effectiveness  of  the  sys- 
2 

tern.  In  any  case,  C_  plays  an  Important  role.  We  now  turn  to  the  quantitative 
2  n 

role  played  by  Cn  in  the  above  mentioned  effects. 

3.2  Amplitude  Fluctuation!  (Scintillation) 

Derivations  of  the  equations  given  below  will  be  found  in  the  references  cited. 

In  the  present  section  we  are  concerned  with  amplitude  variations  of  the  carrier 

signal.  In  the  case  of  these  fluctuations,  a  "log-normal  distribution"  is  involved. 

Tatarski1  (p.  169)  has  shown  that  the  mean  square  fluctuations  of  the  logarithmic 

2 

level  of  the  signal  amplitude,  A,  are  related  to  C 


where  A  is  the  mean  amplitude,  k  is  the  wave  number  of  the  radiation  and  x  is  the 
distance  along  the  path  of  the  beam.  L  is  the  total  length  of  the  beam  path  and  r 
is  the  position  vector.  It  is  assumed  that  /  j  «  Vxi:  «  1°,  that  is,  that  the  quan¬ 
tity  ^ A.I  ,  which  is  known  as  the  size  of  the  first  "Fresnel  zone,"  falls  into  the 
size  range  of  inertial  range  eddies  (I  j  and  i(),  the  reader  will  recall,  are  the 
"inner"  and  "outer"  lengths  respectively). 

It  is  important  to  notice  that  these  scintillation  effects  depend  on  where  they 
occur  along  the  optical  path.  This  is  indicated  by  the  term  in  Eq.  (57)  under 

the  integral  sign,  and  x  is  to  be  interpreted  as  the  distance  from  the  receiver. 

In  the  case  where  «  l  j  (that  is,  the  Fresnel  zone  is  much  smaller  than 

the  inner  scale),  then,  as  Tatarski^  (p-  199)  has  shown 


7.  37  /  1~7/3  J  C2C?)  x2  dx 

o 


(58) 


2 

In  the  simple  case  where  C’n  is  constant  along  the  path,  as  it  might  be  for  the  case 
of  surface-to-surface  communication  for  example,  then  the  latter  may  be  taken  out 
from  tlie  integral  and, 
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(59) 


for  i  j  « 


«  tQ  and 


for  i  j  » 


(60) 


(Ref.  1.  p.  169). 

The  interested  reader  may  wish  to  read  the  next  section  which  gives  a  physical 
2 

explanation  of  Cn  induced  scintillations  in  terms  of  the  random  lens  picture  of 
turbulence.  The  above  formulas  can  be  derived  directly  from  this  simple  picture. 

3.3  Eddie*  a*  Lcnte*  (An  Addendum) 

The  equations  in  subsection  3.  1  can  be  obtained  in  a  very  heuristic  and  visual* 
izable  manner  by  considering  the  turbulent  medium  as  an  array  of  blobs  which  can 
act  as  weak  lenses.  First  we  shall  follow  Clifford^  (pp.  37-39)  and  consider  the 
following  two  questions:  "What  is  the  optimum  size  of  an  eddy,  ).n  relation  to  the 
wavelength  of  radiation,  for  causing  fading  effects?"  and  "How  does  the  weighting 
effect,  x  ;  ,  arise  in  Eq.  (57)?"  After  treating  these  two  questions  we  turn  to  the 
details  of  deriving  the  scintillation  formulae  from  lens-focusing  effects. 

3.3.  1  OPTIMUM  EDDY  SIZE  AND  WEIGHTING  EFFECTS 

Figure  3  depicts  an  eddy  of  dimension  t  perpendicular  to  the  path  of  propaga¬ 
tion  of  an  electromagnetic  signal.  This  path  is  from  left  to  right  in  the  figure  and 
ends  at  the  distance  L.  The  distance  between  the  eddy  (at  z)  and  the  receiver  (at 
L)  is  marked  y  and  is  equal  to  (L  -  z).  For  maximum  destructive  interference 
between  waves  propagating  from  point  P  to  L  and  waves  from  z  to  L,  the  difference 
in  the  path  lengths  P  L  and  y  should  be  one-half  the  wavelength. 

Geometrically 

<PL)2--f2  +  y2  (61) 

and  setting 

(PL  -  y)  «  i  (62) 


11.  Clifford,  S.  F.  (1978)  The  Classical  Theory  of  Wavs  Propagation  in  a  Turbu¬ 
lent  Medium,  in  J.W.  Strohbehn  (editor),  Laser  Beam  Propagation  in  the 
Atmosphere,  Springer  Verlag,  N.Y.  "  . . .  "  ~1“  ~ 
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Figure  3.  Path  Differences  From  Two  Edges 
of  an  Eddy  of  Size  l 


2 

and  eliminating  P  L  in  Eq.  (61)  by  Eq.  (62)  and  ignoring  (X/2)  we  obtain 

<  =  ^X(L  -  z)  (63) 

for  the  first  interference  condition.  Note  that  this  says  t  is  equal  to  "a  Fresnel 
zone  for  the  remaining  length  of  optical  path," 

This  minimum  size,  l ,  for  the  irregularity  is,  it  turns  out,  the  most  effective 
one  for  producing  interference  scintillation.  This  is  true  because  smaller  eddies 
are  less  effective  since  the  turbulent  cascade  is  characterized  by  ever  smaller 
amounts  of  turbulent  intensity  and  fluctuatfon  as  the  wave  number  increases.  As 
Kolmogorov  has  shown  this  dependence  is  k" 3/3  for  the  one -dimensional  spatial  fre¬ 
quency  fluctuation  spectrum.  In  the  case  of  three-dimensional  spectra,  the  index  of 

2 

refraction  fluctuation  spectrum,  $n(k)  depends  on  k  and  Cn  in  the  following  way 
(Ref.  11,  p.  19). 

*n(k)  -  0.033  C2  k"11/3  *  (64) 

(<  j  «  k  *  «  I  )  which  is  a  result  directly  derivable  from 

C*  r2/3  (65) 


*  -5/3 

This  corresponds  to  k  ' 


in  the  one -dimensional  spectrum. 
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From  Eq.  (64),  then,  It  la  clear  that  the  changes  of  index  of  refraction  are  a 
decreasing  function  of  eddy  sizes  and  that,  therefore,  eddies  smaller  than  <  given 
by  Eq.  (63)  are  less  effective  In  causing  effects  on  the  electromagnetic  waves. 

Now  consider  eddies  larger  than  f .  As  eddies  Increase  In  size  another  factor 
enters  which  decreases  their  ability  to  cause  Interference.  As  I  becomes  larger, 
the  blob  becomes  less  effective  in  its  ability  to  cause  interference  because  an 
eddy  of  size  I  can  diffract  a  wave  through  an  angle  of  size  A/i  (in  radians).  How¬ 
ever,  in  order  that  diffracted  light  from  P  reach  L,  it  must  be  diffracted  through 
an  angle  given  by  i  /(L  -  z).  But  since  we  are  considering  an  eddy  larger  than  the 
i  given  by  Eq.  (63),  this  condition  cannot  be  met;  that  Is,  the  larger  eddy  cannot 
diffract  radiation  through  a  large  enough  angle  and  hence  there  Is  no  interference. 
From  these  considerations  we  see  that  I  given  by  Eq.  (63)  causes  the  maximum 
destructive  Interferences  In  the  waves. 

Next,  we  considerW,  the  weighting  function  x®^®  in  Eq.  (57).  The  spatial 
frequency  distribution  for  the  index  of  refraction  inhomogeneities  along  the  optical 
path  is  described  by  the  one -dimensional  Kolmogorov  distribution  which  is,  as 
was  already  mentioned,  dependent  on  the  wave  number  as  k  5^.  Inserting  the 
value  t  =  ^  A(L  -  z)  into  this  form,  that  is,  ,  we  arrive  at  the  weighting 

function  W  given  by 

W~[L-z]5/6  (66) 

This  is  the  physical  reason  for  x5/6  under  the  integral  sign  in  Eq.  (57). 

3.3.2  THE  "WEAK  LENS"  APPROACH  TO  EFFECTS 

12 

In  the  following  we  base  our  treatment  on  that  of  Strohbehn  (p.  72)  but  we 

sacrifice  some  generality  in  the  cause  of  simplicity  and  brevity.  We  now  set  out  to 

derive  Eqs.  (59)  and  (60)  by  simply  regarding  the  blobs  as  simple  lenses  which  can 

be  described  by  geometric  optics.  Figure  4  shows  the  path  length  L,  the  "lens" 

diameter,  i,  and  lens  focal  length  F.  Before  going  any  further  one  might  ask, 

"What  values  of  F  would  these  lenses  have  for  typical  atmospheric  turbulence?" 

12 

According  to  Strohbehn  F  is  in  the  range  of  1000m  to  100,  00C  km.  In  other 
words,  these  turbulent  blobs  are  very  weak  lenses  with  very  long  focal  lengths. 

To  continue,  let  the  refractive  index  of  the  lens  be  given  by  An  =  (n  -  1).  For  such 
a  lens,  with  diameter  l,  the  focal  length  can  be  taken  as 

12.  Strohbehn,  J.W.  (1978)  Modern  theories  in  the  propagation  of  optical  waves 

inJ.W.  Strohbehn  (editor),  Laser  Beam  Propagation  in  the  Atmosphere. 
Springer  Verlag,  N.  Y. 

13.  Jenkins  and  White  (1957)  Fundamental  of  Optics  (3rd  Ed.  ),  McGraw-Hill, 

N.  Y. 
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Figure  4.  Eddy  Lens  and  Scintillation 


F  = 


i 

An 


(67) 


(cf.  Ref.  13,  p.  51).  When  An  is  positive,  the  lens  wili  focus  the  radiation  and 
when  negative  it  will  defocus  it.  We  now  use  geometrical  optics  to  calculate  the 
variations  in  the  amplitude  due  to  these  lenses.  To  do  this  we  make  use  of  con¬ 
servation  of  energy  in  the  beam.  In  Figure  4  we  consider  the  case  where  the  blob 
lens  of  size  i  focuses  a  beam  of  width  f  down  to  a  width  2r  at  point  R  (location  of 

receiver).  Assuming  that  the  original  amplitude  of  the  beam  is  A  ,  and  using  the 

2  ° 

fact  that  the  radiant  energy  goes  as  A  (per  unit  of  area),  conservation  of  energy 
leads  to 


( 

I 

l 

I 

i 

! 

i 


aj  nr  -  a2  -z 


(68) 


where  A  is  the  amplitude  at  point  R,  Now  we  define  6A  s  (A  •  Aq).  Solving 
Eq.  (68)  for  6 A  we  have 


‘A-ao  (£  -  l) 


(69) 


From  Figure  4  we  can  see  that  the  triangle  with  height  i/2  and  base  F  has  the 
same  angle  (at  the  extreme  right  of  the  figure)  as  the  triangle  with  height  r  and 
base  (F  -  (L  -  z)| .  The  tangent  of  this  angle  gives  us 


r 

F  -  (L  -  z) 


(70) 
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Solving  this  for  l/2r  anil  Inserting  the  result  In  Eq.  (GO)  we  arrive  at 


**■ 

o 


F  -  (L  -  7) 


(71) 


where  we  used  (L  -  z)  «  F  as  an  approximation. 

Next,  we  consider  an  array  of  such  lenses  with  Indices  of  refraction  which 
depend  upon  the  sizes  In  the  manner  of  the  Inertial  range  structure  function 
(Eq.  (42)) 


<An 2(f)>  ~  C2  I2/3 


(72) 


where  we  have  Identified  O'2  with  (An2)  and  C2  with  C2.  Of  course,  Eq.  (72)  is 
correct  only  when  I .  «e  I  «  f  .  In  the  case  where  t  «  t ,,  It  can  be  shown  that 

1  o  n  O  1 

one  must  replace  i  by  t  /t  J  In  Eq.  (72).  Using  Eq.  (72)  In  Eq.  (67) 


F 


(73) 


where  we  define 


i '  ■  i  ,  (f  j  «  i  «  iq) 

*'-4  .  <* «*i> 

*i 


(74) 


From  Eq.  (71)  then, 


«') 


nl/3 


(L  -  z) 


Thus 


I4'  C„VZ/3<L-’>  ■  «  <<!» 

o 

Cn  f  "2//3  (L  -  z)  ,  (f  J  «  f  «  /q) 


(75) 


(76) 
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If  we  titume  that  i  ^  ta  larger  than  the  value  t  ■  >J(\L)  (the  most  effective 
scattering  length),  then  In  general  we  can  make  the  approximation  that  eddies  of 
size  1 1(  do  all  of  the  scattering  or  focusing  etc,  and 


C2  t 
<~n  *  1 


-"/3  (L  . 


(771 


o 

Next  we  must  calculate  the  average  value  of  ( 6A/Aq )  .  The  number  of  Inhomo- 
genelttes  along  the  optical  path  we  shall  assume  is  N  -  L/f  In  effect  we  are 
using  a  model  where  all  the  eddies  are  of  size  i  1  and  touching  each  other.  Then 

77  " 


2  2 
where  (L  -  z)  is  the  average  over  the  path.  Taking  this  as  L  we  arrive  at 


6A 

1 r 


C2  L3i  -?/3 
n  l 


(79) 


But  since 


jp  =  6(fn  A)  *  in  A  -  in  Aq 


we  can  write  Eq.  (79)  as 


i 


-7/3 

i 


where 


(t .  «  I  «  I  )  and  (i  <  t  ,) 

t  o  c.  1 


(80) 


(81) 


This  is  Eq.  (60)  which  has  been  derived  (originally  by  means  of  rigorous  analytic 
methods).  The  above  simple  model  arrives  at  the  same  result  with  all  but  the 
constant  of  order  unity. 

In  the  alternative  case  of  1  c  >  /  ^  one  must  replace  /  ^  by  yj  XL  since  in  this 
case  t  is  the  most  efficient  size  of  eddy,  and 
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(82) 


< 


c2  k?/«  L“/« 


where 

*c  >  t  j  (or  1 1  <  V<XL))  . 


This  is  the  same  as  Eq.  (59)  and  again  we  have  everything  but  the  numerical  factor 
from  the  simple  model. 

3.3.3  SATURATION  EFFECTS 

As  more  eddies  enter  the  beam,  multiple  scattering  can  occur  whereby  a 

scattered  ray  can  be  rescattered.  As  this  process  becomes  more  dominant,  the 

phase  of  the  radiation  across  the  beam  becomes  more  and  more  Incoherent.  When 

the  point  is  reached  where  there  is  no  longer  any  coherence  across  a  distance  of 

the  size  of  the  first  Fresnel  zone  (^ XL)  then  the  above  treatment  breaks  down, 

and,  as  a  result,  the  scintillation  effects  become  no  longer  dependent  on  an  increase 

of  path  length.  This  phenomenon  is  known  as  "saturation"  and  it  has  received  a  lot 

of  attention  in  the  literature.  A  detailed  explanation  of  this  in  terms  of  random 

12 

lenses  will  be  found  in  Strohbehn  p.  79. 

3.3.4  THE  DIFFRACTION  GRATING  PICTURE 

For  completeness  it  should  be  mentioned  that,  in  addition  to  the  "blob  lens" 

o 

approach,  there  is  another  simple  picture  available  for  describing  C  effects  due 

14  n 

to  Tatarski,  p.  114.  This  is  a  model  based  on  the  idea  that  a  random  array  of 

diffraction  gratings  of  random  spatial  frequencies  can  explain  the  formula  for  the 

"scattering  cross-section"  of  turbulent  volumes  of  atmosphere.  The  eddies  have 

not  only  random  periods  but  random  orientations.  The  "grating"  picture  arises 

from  the  fact  that  this  particular  theory  makes  intrinsic  use  of  the  spatial  frequency 

spectra. 


3.4  Pliaae  and  Ray  Angle  Fluctuations 

As  first  shown  by  Tatarski/1  p.  170,  and  mentioned  in  Pratt,  111  p.  135,  the  var¬ 
iance  in  phase  as  a  function  of  separation,  p,  across  a  beam  cross-section,  that  is 
to  say  the  structure  function  of  the  phase  ( <j> )  fluctuations  between  two  points  of 
separation  p,  i3  given  by 


14.  Tatarski,  V.  I.  (1971)  The  Effects  of  the  Turbulent  Atmosphere  on  Wave 
Propagation,  National  Science  Foundation,  TT-68-50464. 
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ragiW.ii.iJB 


^(p)  = 


1.46 

2.91 


L 


5/3  f  cl{r)  dx  1 

o 

(i  J  «  P  <K  ^  AL  ) 

L 

(83) 

o 

(t  o»  p  »  yj  AL ) 

The  integration  here  is  carried  out  along  the  ray  starting  from  the  receiver  or 

2 

observing  point  to  the  source.  In  the  case  where  C  is  constant  along  the  path,  the 

2  ^ 

Integrals  can  be  replaced  by  C  L. 

Phase  fluctuations  are  caused  by  two  physically  different  processes.  In  one 
case  there  are  the  effects  of  wavefront  distortions  over  the  beam  cross-section. 

In  the  other  case  there  is  the  effect  of  beam  steering  or  image  dancing  where  the 
wavefront  is  tilted  through  the  angle  (3,  where 


(84) 


where  <j>  is  the  phase  of  the  radiation  (Ref.  10,  p.  135). 

3.5  Spatial  Coherence  Degradation 

One  of  the  important  parameters  in  the  design  of  optical  receivers  is  the 

parameter  rQ  (Pratt,  p,  140-141  and  Fried1*’).  This  quantity,  rQ,  is  a  sort  of 

coherence  length.  Physically  rQ  is  the  distance  such  that,  if  an  antenna  diameter 

is  increased  beyond  rQ,  there  is  significantly  less  improvement  in  performance. 

In  other  words,  rQ  represents  the  coherence  across  the  beam  cross-section.  As 

the  turbulence  degrades  the  phase  of  the  radiation,  r  is  reduced  in  size.  As  C2 

15  °  n 

becomes  larger,  rQ  becomes  smaller.  As  Fried  has  shown  (by  means  of  the 
structure  functions  of  both  amplitude  and  phase  fluctuations) 

rQ  =  [i.  2  X  10"8]  (A)6/5  L'3/5  C"6/5  (85) 

where  L  is  the  path  length,  constant  Cn  is  assumed  along  the  path,  and  A,  as  usual, 
is  the  wavelength  of  the  radiation. 

In  addition  to  the  determination  of  the  optimum  aperture  size,  the  value  of  r 

1  fi  O 

is  used  to  calculate  the  resolution  of  focused  images.  It  is  also  an  important 

15.  Fried,  D.  L.  (1967)  Optical  heterodyne  detection  of  an  atmospheric  distorted 

signal  wave  front,  Proc,  IEEE  55:57. 

16.  Beran,  M.J.,  and  Parrent,  G.  B.  (1964)  Theory  of  Partial  Coherence. 

Prentice -Hall,  N.J. 


parameter  In  the  calculation  of  so-called  "saturation  effects"  (that  is,  where  an 
increase  of  the  path  length  of  the  ray  no  longer  increases  the  value  of 
«4A/Ao)2}). 

From  a  practical  point  of  view  in  the  context  of  systems  design  (where  adaptive 
optics  is  employed)  the  most  Important  application  of  rQ  Is  to  the  determination  of 
the  number  of  activator  elements  needed  to  achieve  a  given  performance.  Thus  rQ 
determines  the  number  of  segments,  for  example,  in  a  compensating  mirror.  If 
rQ  is  too  small,  it  might  imply  that  the  present  state  of  the  art  of  electronic  optical 
compensation  is  insufficient  for  certain  applications)  thus,  this  parameter  is  in¬ 
deed  of  paramount  importance. 

It  is  of  interest  to  understand  the  physical  significance  of  r  .  Figure  5  shows 
this  dlagrammatically.  Consider  two  rays,  1  and  2,  which  travel  from  left  to 
right  in  the  figure.  In  Figure  5A,  these  rays  are  separated  by  a  distance  r  <  rQ  in 
the  plane  of  the  beam  cross-section.  The  value  of  rQ  corresponds  to  the  size  of 
the  "eddy  lens"  affecting  the  phases  of  the  ►•ays.  As  can  be  seen,  since  both  rays 
go  through  the  same  "lens,  "  one  can  assume  (in  this  model)  that  they  both  suffer 
the  same  change  in  phase  and  hence  maintain  their  phase  relation  after  the  inter¬ 
action.  The  alternative  case  is  shown  in  Figure  5B.  Here  r  >  r  and  the  two  rays 
go  through  different  lenses  of  different  index  of  refraction.  This  would  destroy 
their  original  phase  relation.  After  a  number  of  such  encounters  the  two  rays 
would  be  "incoherent"  with  respect  to  each  other.  This,  in  essence,  is  the  phys¬ 
ical  meaning  of  rQ.  The  following  description,  which  is  based  on  a  number  of 
papers  by  Yura,17  will  render  the  above  qualitative  description  more  quantitative 
and  at  the  same  time  explain  Eq.  (85). 

Let  ASj  be  the  change  in  phase  between  rays  1  and  2  (initially  with  the  same 
phase)  caused  by  going  through  an  inhomogeneity  where  r,  again,  is  the  distance 
between  the  rays.  Let  An  be  the  difference  of  index  of  refraction  between  the  two 
eddies  traversed  by  the  rays.  ■  Then  the  difference  in  phase,  AS^  caused  by  this 
interaction  is 

ASj  =  kAn  •  r  (86) 

In  Eq.  (86)  it  is  assumed  that  r  is  also  the  diameter  of  each  of  two  round  eddies  as 
in  Figure  CB.  Note  that  r  An  is  the  optical  path  difference  and  that  the  subscript, 

1,  denotes  a  single  interaction. 

If  one  recalls  the  definition  of  the  structure  function  (Eq.  (2))  one  obtains 

(An2)  =  ^(r)  (87) 


17.  Yura,  H.T.  (1974)  J.  Opt.  Soc.  Am.  64:59,  1526,  357,  1211. 
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A 


Eddy.  B:  Two  Rays  Through  Two  Eddies 


(r  is  the  diameter  of  the  "lens"  or  ~  outer  length)  thus 

(AS2)  -  k2  r2  ^(r)  .  (88) 

^  2  2  /3 

But  in  the  inertial  range  we  have  -2r(r)  ~  C*  r  '  ,  thus 

(AS2)  -  k2  r8/3  C2  .  (89) 

We  must  now  consider  multiple  scattering.  The  number  of  encounters,  N, 
with  "lenses,  "  along  a  path  of  length  L  (between  source  and  receiver)  can  be 
estimated  by 

N  =  h  (90) 

r 

2 

where  we  have  assumed  homogeneity  for  simplicity.  In  order  to  obtain  (AS  (r)) 

due  to  N  events,  we  must  remember  that  the  phase  changes  would  occur  in  a  ran- 

2  2 
dom  manner.  For  this  reason  (AS  (r))  would  increase  as  N  rather  than  N 

(Jenkins  and  White,  *3  p.  218).  Thus 

(AS2(r)>  =  N(AS2(r)>  ~  k2  C2  L  r5/3  .  (91) 
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When  r  =  r  we  assume  that  there  Is  incoherence.  For  the  latter  we  take 
o 

<AS2(r0)>  =  n  2  (92) 

and  insert  this  into  Eq.  (91)  to  obtain 


,2  . 


3/5 


(93) 


which  is  within  a  constant  of  Eq.  (85). 

Finally  we  conclude  with  a  very  simple  description  of  why,  physically,  rQ 
determines  the  maximum  aperture  mentioned  previously  in  connection  with  image 
resolution.  This  explanation  is  due  to  George  Parrent  (private  communication). 
First  it  must  be  recalled  that  the  resolution  of  an  image  is  determined  by  the  dif¬ 
fraction  pattern  of  the  effective  aperture.  The  latter  is  controlled  by  the  destruc¬ 
tive  interference  between  the  two  most  distant  parts  of  the  incident  wavefront.  But 
interference  (or  the  concept  of  wavefront  for  that  matter)  demands  the  presence  of 

coherence.  As  we  have  seen  r  is  the  coherent  limit;  therefore  r  is  the  effective 

o  o 

size  of  the  aperture  so  far  a3  interference  is  concerned.  If,  for  example,  an 
aperture  had  a  radius  of  4  rQ,  it  would  gather  much  more  light  than  one  of  radius 
r  ;,  regarding  interference  effects  it  is  as  if  there  were  several  different  lenses 
operating  with  different  wavefronts  without  any  mutual  interference  effects  so  that 
the  diffraction  pattern  would  be  that  of  one  lens  of  radius  r  .  It  is  now  obvious 
why  an  aperture  larger  than  rQ  would  not  help  much  with  resolution  problems. 


3.6  Beam  Spreading 


Since  a  light  beam,  or  any  electromagnetic  beam,  does  not  have  sharp  edges 

(due  to  diffraction)  its  width  is  often  defined  by  its  Gaussian  (1/e)  radius.  Desig- 

1 8 

nating  this  width  by  agj.j  it  has  been  shown  (Ishimaru,  p.  137)  that 


2 

aeff 


a2  +  1  6  C2  k*/2 
o  n 


1/3  8/3 
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aeff  " 


tUC2^3 
n  o 


(n  »  i)  i 

(n  «  i)  ! 
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18.  Ishimaru,  A.  (1978)  The  beam  wave  case  and  remote  sensing  in  J.W. 

Strohbehn  (editor).  Laser  Beam  Propagation  in  the  Atmosphere,  Springer 
Verlag,  N.  Y. 
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where  L  Is  path  length,  aQ  is  Initial  Gaussian  width  and  n  ■  (aL)"1  where  a  is 
defined  by 


a  ■ 


(95) 


2  2 

Notice  that  as  CR  increases,  increases  linearly. 

3.7  Typical  Values  -  Optical  Compensation  Techniques 

2 

In  the  above  subsections  we  have  seen  quantitatively  how  CR  affects  design 

parameters  of  optical  systems.  We  now  complete  Section  3  by  considering  two 

2 

additional  related  topics,  namely:  (a)  what  are  typical  values  of  Cn  in  the  atmos¬ 
phere,  and  (b)  what  are  some  of  the  techniques  available  for  the  purpose  of  mini- 

2  19 

mizing  Cn  effects  on  systems.  According  to  Davis 


Weak  turbulence 

Cn  =  8  X  10'9  m*1/3* 

(96) 

Intermediate  turbulence 

C  =  4  X  10"8  m-1/3 
n 

(97) 

Strong  turbulence 

C  =  5  X  10'7  m'1/3 
n 

(98) 

These  are  meant  to  characterize  average  daytime  conditions  a  few  meters  above 

the  ground.  C  decreases  roughly  as  the  -1/3  power  of  altitude  within  the  first 
n  19 

100  m  of  the  surface.  Hufnagel  and  Stanley  are  quoted  by  Davis  as  having  inte- 
o 

grated  Cn  along  a  vertical  path  through  the  atmosphere  under  daytime  turbulence 
conditions.  They  found 

C  -  1.3  X  10'11  m’1/3  .  (99) 

Clifford*1  has  quoted  the  following  values  from  Obukov.  In  the  daytime, 

84  percent  of  observations  fell  within 

2.  3  X  10'7  m"1/3  <  C  <  7.  4  X  10'7  m"1/3  .  (100) 

n 

At  night,  59  percent  (note,  a  larger  spread  due  to  the  smaller  percent)  fell  within 
7.4  X  10'8  m"1/3  <  Cn  <  2.3  X  10'7  m'1/3  (101) 


*m  for  meters. 

19.  Davis,  J.  I.  (1967)  Consideration  of  atmospheric  turbulence  in  laser  systems 
design,  Applied  Optics  5:139-141). 
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Clifford**  also  quotes  Neff's  published  values 


(C2)1^2  ~  2.  f>  X  10'8  m-1,/3  ±  3.  9  X  10'8  m"*/3 


(102) 


for  early  morning,  and 

,2-  */2 


<c‘> 


5.  5  X  10~8  m'*/3  ±  1.4  X  10“7  m‘*/'3 


(103) 


for  afternoon  measurements. 

20 

Shapiro  has  described  some  of  the  compensation  techniques  used  to  minimize 
o 

C  problems.  The  following  are  some  brief  descriptions  of  these. 

Interferometry:  Consider  the  classical  problem  of  measuring  the  distance  be¬ 
tween  two  very  distant  stars.  Such  a  measurement  would  be  impossible  on  the 

2 

earth's  surface,  due  to  C  effects,  if  tiie  latter  were  not  somehow  eliminated.  It 

n 

turns  out,  however,  that  the  image  of  the  interference  pattern  caused  by  two  slits 
is  not  displaced  by  turbulence-induced  phase  distortions.  An  interference  pattern 
is  the  autocorrelation  of  the  optical  signal.  From  this  fact  one  can  derive  the  re¬ 
sult  that  it  is  phase  independent.  The  general  problem  of  measuring  the  distance 
between  two  seif-luminous  objects  can  be  solved  in  this  manner  up  to  limits  which 
do  not  depend  on  turbulence. 

Phase  compensated  imaging:  As  surprising  as  it  may  seem  to  the  uninit  iated, 

it  is  now  possible  to  actively  measure  and  compensate  for  turbulence-induced  phase 

perturbations  in  real  time  (to  some  useful  degree).  This  is  made  possible  by  the 

21 

fact  that  atmospheric  phase  modulations  occur  at  a  rate  less  than  1  kHz,  and 
that  recent  advances  in  optical  and  electronic  technology  can  enable  the  compensa¬ 
tion  to  be  made  automatically.  The  device  which  does  this  is  made  up  of  a  wave 
sensor  (or  phase-estimator  array),  a  wavefront  corrector  (deformable  mirror) 
and  a  feedback  system  which  adjusts  the  corrector  in  a  manner  which  maintains  e 

constant  image.  The  term  "rubber  mirror"  has  been  used  to  describe  this  device 

■’0 

and  the  reader  will  find  a  number  of  references  in  Shapiro. 

20 

Diversity  reception:  Fading  can  be  minimized  by  using  a  number  of  receiv¬ 

ers  which  are  sufficiently  separated  (so  that  they  occupy  different  parts  of  the 
reception  pattern)  and  by  combining  their  signals  in  an  optimal  manner.  This  is 
used  in  downlink  space  to  earth  channels  and  is  a  well  known  technique  in  commun¬ 
ication  engineering.  For  obvious  reasons  it  cannot  be  used  to  dec  reuse  scintillation 
effects  in  the  uplink  channel. 

20.  Shapiro,  J.H.  (1978)  Imaging  and  optical  communication  through  atmospheric 

turbulence,  inJ.VV.  Strohbehn  (editor),  baser  Beam  Propagation  in  the 
Atmosphere,  Springer  Verlag,  N.Y,  -----  - 

21.  Greenwood,  D.  P.  (1977)  Bandwidth  specification  for  adaptive  optics  systems. 

Optical  Soc.  of  Am.  87:390-392. 
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Adaptive  variable  rate  communication:  Kennedy  et  al  have  suggested  an 
idea  that  may  have  use.  As  has  been  already  pointed  out,  the  amplitude  fluctua- 
tions  are  log-normally  distributed.  This  implies  that  the  distribution  is  highly 
skewed  and  it  has  been  observed  that  amplitude  fluctuations  that  far  exceed  the 
mean  are  not  improbable.  For  example,  astronomers  report  that  with  small 
apertures  sudden  increases  of  intensity  of  400  percent  above  mean  values  over  a 
time  of  1/100  sec.  are  not  unusual.  By  adjusting  the  bit  rate  of  transmission  it 
may  be  possible  to  take  advantage  of  such  periods  of  high  gain.  The  practicality 
of  such  a  scheme  seems  not  to  have  been  investigated,  however. 


4.  STATE  OF  THE  ART  FOR  ESTIMATION  OF  df  FROM 
RAWINSONDE  DATA 

4.1  State  of  the  Art 

3 

VanZandt  et  al  have  devised  a  method  to  relate  standard  raw  insonde  measure  - 

o 

rnents  to  direct  radar  measurements  of  C  .  If  this  leads  to  a  comprehensive  way 
2  ^ 

to  deduce  Cn  (for  most  altitudes  of  interest)  from  already  existing  rawinsonde  data, 
it  will  be  possible  to  use  the  large  data  bank  of  the  latter  in  order  to  ascertain 

9 

geographical,  seasonal,  diurnal  and  synoptic  weather  pattern  effects  on  Cn.  As  we 

shall  see,  however,  some  problems  remain  to  be  solved  before  such  is  possible. 

3 

The  work  of  VanZandt  et  al  is  the  second  version  of  the  model  that  they  have 
proposed.  The  theory  behind  this  model  is  based  in  part  on  a  model  of  clear  air 
turbulence  which  assumes  that  the  shear  instability  is  the  sole  cause  of  the  turbu- 
lence.  Some  of  the  work  by  Rosenberg  and  Dewan  was  incorporated  into  their 
model.  While  this  latter  paper  was  originally  intended  primarily  for  describing 
turbulence  in  the  stratosphere  it  is  also  appropriate  for  the  modeling  of  clear  air 
turbulence  in  the  troposphere  (known  popularly  as  CAT).  The  main  difference 
between  CAT  in  the  troposphere  and  stratospheric  turbulence  is  that,  in  the  trop¬ 
osphere,  the  stability  is  significantly  lower  (that  is,  the  buoyancy  frequency  is 
down  by  about  a  factor  of  two).  In  addition,  there  are  other  causes  of  turbulence 
in  the  troposphere  having  to  do  with  convection  and  weather  fronts,  etc.  At  the 
lowest  altitudes  (of  order  1  km)  there  is  the  so-called  boundary  layer  where  tur¬ 
bulence  effects  depend  on  factors  relating  to  surface  influences. 

3 

In  order  to  explain  the  model  of  VanZandt  et  al  it  is  necessary  to  review  the 
nature  of  stratified  turbulence.  (A  much  more  comprehensive  treatment  will  be 

22.  Kennedy,  R.S.,  and  Karp,  S.  ( 1969)  Optical  Space  Communication,  NASA 

SP-217.  - - 

23.  Rosenberg,  N.W.,  and  Dewan,  E.M.  (1975)  Stratospheric  Turbulence  and 

Vertical  Effective  Diffusion  Coefficients.  AFtRL-TR -75-0519, 

AD  A019  708,  Air  Force  Geophyoics  Lab.,  Bedford,  MA, 

♦This  may  be  disputed  somewhat,  but  it  does  not  affect  the  argument. 
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found  In  Rosenberg  and  Dewan.  )  The  mechanism  for  turbulent  breakdown  in 
stratified  fluids  Is  the  Kelvin-Helmholtz  Instability  (KH).  The  crucial  parameter 
In  the  theory  of  this  instability  is  the  Richardson  number  defined  by 

o  g  dF/dz  g  (dT/dz  +  y) 

n-  -  — -  — n - ty  -  Tpr-  - - - «-  (104) 

1  7  (dU/dzr  %  <dU/dz)2 

where  0  is  the  potential  temperature,  T  is  the  absolute  temperature,  y  is  the  adia¬ 
batic  lapse  rate,  U  is  the  average  horizontal  velocity  and  the  subscript  naught 
indicates  a  larger  scale  average  over  the  region  of  interest.  When  R,  <  1/4  It  Is 
possible  (but  not  mathematically  guaranteed)  that  turbulence  breakdown  will  occur; 
in  other  words,  it  has  been  shown  that  R;  <  1/4  ts  a  necessary  but  not  sufficient 
condition  for  turbulence.  In  nature  it  seems  safe  to  assume  that  it  is  generally 

also  a  sufficient  condition.  25,  26  For  a  simple  energetic  arguement  for  why  this  is 

27 

true,  the  reader  may  consult  Businger  who  demonstrated  that  R^  =  1/4  is  the 
condition  where  the  kinetic  energy  (available  from  the  shear  across  a  layer  of 
given  thickness)  is  equal  to  the  work  necessary  to  overcome  the  buoyancy  force 
(due  to  atmospheric  stability)  in  the  process  of  overturning  the  lEyer.  Figure  6 
illustrates  the  sequence  of  configurations  (from  the  generation  of  a  wave  to  the 
turbulent  breakdown)  in  the  development  of  the  K-H  instability.  Notice  that  before 
the  breakdown  occurs,  the  wave  "rolls  up."  The  actual  final  breakdown  seems  to 
be  due  to  the  fact  that  when  roll  up  occurs  there  are  layers  of  dense  fluid  over  less 
dense  fluid  and  numerous  local  convective  "turnover"  events  simultaneously  occur 
within  these  superposed  layers. 

28 

Figure  7  shows  a  diagram  baaed  on  the  one  given  by  Browning  and  Watkins. 
Those  investigators  worked  with  radar  observations  of  CAT  in  the  lower  tropo¬ 
sphere.  It  may  be  of  interest  that,  such  radar  sightings  were  historically  known 


24.  Miles,  J.W.  (1961)  On  the  stability  of  heterogeneous  shear  flows,  parti, 

J.  Fluid  Mech.  10:496-512. 

25.  Woods,  J.D.  (1968)  An  investigation  of  some  physical  processes  associated 

with  the  vertical  flow  of  heat  through  the  upper  ocean.  The  Meteorological 
Mag.  97:65-72. 

26.  Thorpe,  S.A.  (1973)  Turbulence  in  stably  stratified  fluids:  A  review  of 

laboratory  experiments,  Boundary  Layer  Met.  5:94-119. 

27.  Businger,  J.A.  (1969)  On  the  energy  supply  of  clear  air  turbulence,  in 

Y.  H.  Pao  and  A,  Boldburg  (editors).  Clear  Air  Turbulence  and  its 
Detection,  Plenum  Press,  100-108. 

28.  Browning,  K.  A.,  and  Watkins,  C.  D.  (1970)  Observations  of  clear  air  turbu¬ 

lence  by  high  power  radar.  Nature  227:260-263. 


Figure  6.  Development  of  Billow  and  Subsequent  Turbulence 
(After  Woods) 
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Figure  7.  Browning  and  Watkins'  Description  of  Radar 
History  of  CAT 
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as  "braid  angels"  before  their  nature  was  understood.  Atlas  et  al  seem  to  be 
the  first  to  have  identified  these  as  CAT. 

Before  the  turbulence  commences  there  i3  a  difference  of  potential  tempera¬ 
ture  across  the  layer.  This  is  caused  by  the  fact  that  the  atmosphere  (In  the 
absence  of  convection  caused  by  heating  from  below)  is  usually  stable.  When  tur¬ 
bulence  takes  place,  the  mixing  causes  parcels  of  air  that  were  Initially  at  differ¬ 
ent  altitudes  (having  different  potential  temperatures)  to  be  moved  together.  This 
causes  the  index  of  refraction  to  vary  wildly  (in  comparison  to  the  original  smooth 

configuration)  and  hence  to  scatter  the  radar  signal  back  to  the  receiver.  In  this 

2 

manner  radar  can  measure  the  C  in  its  beam. 

2  n 

To  relate  Cn  to  raw  insonde  data  one,  ideally,  would  plot  R.  qs  a  function  of 
altitude  and  infer  from  such  stability  profiles  an  estimate  of  the  turbulent  layer 
configurations  (that  is,  by  designating  as  turbulent,  or  about  to  be  turbulent,  those 
regions  where  Rj  <  1/4). 

Of  course,  the  rawlnsonde  measurements  have  low  resolution  —  too  lew  In  fact 
to  directly  measure  the  turbulence  configuration  directly. 

29.  Atlas,  D. ,  Metcalf,  J.I.,  and  Gossard,  E.  E.  (1970)  The  birth  of  "CAT"  and 
microscale  turbulence,  J.  Atm.  Sci.  27:903-913. 


To  actually  obtain  an  estimate  of  C2,  one  makes  use  of  Eq.  (54)  repeated 
here  for  convenience. 

Cj  •  (CONST.)  IQ4/3  M2  (54) 

2 

where  M  can  be  estimated  from  the  rawlnsonde  data  (which  gives  temperature 
and  specific  humidity  as  a  function  of  altitude).  Notice  in  Eq.  (54)  there  Is  an 
unknown,  i^.  The  constant  Is  known  provided  that  the  ratio  of  eddy  dlffusivttles 
Is  taken  equal  to  unity. 

The  radar  estimate  for  C2  Is  obtained  from3 

n 

C2(radar)  «  19.  7  X  10“30  [Tc  +  1875]  r2  (  |  )  [m*2/3]  (105) 

where  r  is  range,  (S/N)  is  the  signal -to-nolse  ratio,  and  Tc  Is  the  cosmic  noise 
temperature.  By  comparing  results  from  Eq.  (105)  and  Eq,  (54)  they  were  able 
to  refine  Eq.  (54)  in  a  manner  which  made  it  useful;  but,  first  they  had  to  solve 
certain  problems. 

The  radar  observations  had  a  resolution  of  about  1  km.  In  contrast  one  would 
expect  the  turbulent  layers  to  have  an  Individual  thickness  of  one-tenth  this  size. 
Another  problem  is  introduced  by  the  fact  that  the  wind  and  temperature  are  tele¬ 
metered  at  intervals  of  about  100  m.  Furthermore  it  was  found  necessary  to 
average  this  data  over  several  such  altitudes  in  order  to  obtain  reliable  data.  By 
means  of  such  running  averages,  profiles  of  average  shear,  (S*),  were  obtained. 

As  already  mentioned  such  shears  would  not  reveal  any  regions  where  R^  <  1/4. 

To  overcome  this  problem,  VanZandt  et  al3  inferred  the  existence  of  "micro- 
shears"  or  "fine  structure  shears"  which  would  cause  turbulent  layers  to  form. 

They  devised  a  formalism  to  relate  the  population  of  microshears  to  (S) .  In  order 
to  do  this  they  used  the  following  simple  and  ingenious  technique.  First  they 
assumed  that  the  values  of  the  microshears  were  normally  distributed  about  (S) 
and  that  the  variance  of  this  distribution  was  a  fixed  quantity.  Thus,  starting  with 
a  measured  (S)  profile  it  was  possible  for  them  to  infer  the  statistical  probability  of 
the  occurrence  of  turbulence  layers.  In  particular,  they  could  infer  the  fraction  of 
altitude,  occupied  by  a  given  slab,  which  is  expected  to  be  in  the  state  of  turbulence 
on  the  basis  of  <S)  and  the  value  or  the  standard  deviation,  cr.  Like  l  above,  a  was 
initially  an  unknown  quantity.  Since  the  values  of  <S>  were  usually  far  too  small  to 
cause  turbulence,  it  was  only  the  shears  in  the  "tails"  of  the  Gausstan  distributions 
which  would  cause  R.  <1/4  and  hence  turbulence.  As  will  be  discussed  below, 


Caution;  this  symbol  wao  previously  used  for  phase. 
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V.i.  .  -  ruw.u.i., 


wn 


a  and  t  were  determined  by  means  of  an  optimal  fit  with  the  data  when  rawlnsonde 

0  2 
were  compared  to  radar  measurements  of  C  . 

In  order  to  make  possible  the  comparison  with  1-km  resolution  radar  meas- 

2 

urements,  they  calculated  the  mean  values  of  CR  within  the  slabs,  labeled  "l",  by 
means  of 


C 


n  t 


Cn,  l  F1 


(106) 


2  2  2 

where  C  ,  is  the  mean  value  of  C  In  turbulent  layers  In  slab  1  and  C  .  Is  the 
n,  i  n  n,  l  „ 

mean  value  Including  both  the  turbulent  layers  and  non-turbulent  regions.  CT  t  is 

determined  from  rawlnsonde  data  by  Eq.  (54)  and  by  Eq.  (105),  and  F.  Is  the  mean 

2  1  2 

fraction  of  the  slab  that  Is  turbulent.  Calling  (C‘  (model))  the  mean  square  of  C* 

over  a  kilometer  (In  order  to  make  a  comparison  with  1-km  resolution  radar),  we 

have 


<c;(mod.t»  ■  2  <  i  [fi^i]  <>«> 

i 

t  h 

where  Az^  is  the  thickness  of  the  i  n  slab  and  where  the  sum  is  taken  so  that  the 
Az.  adds  to  1  kilometer.  This  leads  to 

(Cjlmodel))  .  £  a2  lQ4/3  M.  F.  [~ ^  (108) 

i 

2 

where  a  is  about  2.  8  and  a'  1.  Both  the  value  of  I  and  a  are  determined  (once 

o 

and  for  all)  by  an  optimum  fitting  to  the  data  which  involves  a  careful  comparison 
2  2 

between  (Cn  (model))  with  (Cn  (radar)),  the  latter  being  obtained  via  Eq.  (105)  and 
measurements  of  (S/N)  from  radar  returns.  The  radiosonde  and  radar  measure¬ 
ments  are  made  as  close  together  in  space  and  time  as  is  practical  in  order  to 

3 

make  the  appropriate  estimates  of  t  and  a,  and  VanZandt  et  al  found  l  =  10  m, 

o’  ,  o  ’ 

and  a  =  0.010  s  L  for  the  troposphere  and  t  -  10  m,  a  =  0.015  s  1  for  the  strato¬ 
sphere.  These  values  of  o  are  almost  exactly  the  same  as  those  of  Rosenberg 
23  - 1 

and  Dewan,  for  example,  ct  =  0.014  s  was  given  there  for  the  stratospheric 
case.  As  for  their  value  of  i  =  10  m,  this  "outer  length"  was  close  to  that  used 
by  Crane  on  the  basis  of  isotropy  (that  is,  i  for  Crane,  was  its  value  when  the 
turbulence  became  anisotropic:  private  communication).  The  value  of  I  =  10  m 

*  th 

M  is  defined  by  Eq.  (56)  and  subscript  "i"  refers  to  i  layer. 


43 


Is  of  order  1  /10th  of  the  thickness  of  sn  average  turbulent  layer  in  the  stratosphere 

30 

(aa  seen  by  balloon  measurements.  Barat,  and  by  R,  number  profiles,  Rosenberg 

23  1 

andDewan,  )  which  is  about  200  m.  Indeed  one  would  expect  10  m  to  be  near  the 

31  3 

scale  where  isotropy  breaks  down  ;  but,  In  the  case  or  VanZandt  ot  al,  the  value 

2 

of  I  may  actually  have  more  to  do  with  the  scale  at  which  high  gradients  of  CR 
would  be  expected  to  occur. 

The  main  question  that  now  arises  Is:  "How  well  does  the  model  of 

VanZandt  et  al  actually  work?"  In  their  own  words,  "The  resulting  profiles  of 

2  2 
(C‘  (model))  agree  well  with  the  measured  profiles  of  (CT  (radar))  in  general  shape, 

in  changes  from  day  to  day,  and  in  many  details  from  km  to  km.  This  agreement 

2 

Implies  that:  (1)  The  vertical  profile  of  C  can  be  measured  by  Doppler  Radar, 

2  “ 

(2)  The  vertical  profile  of  Cn  can  also  be  estimated  by  calculation  from  routine 
rawtnsonde  profiles,  using  our  theoretical  model."  They  also  pointed  out  that  the 
agreement  in  absolute  values,  on  the  other  hand  (while  very  good)  only  implies 
that,  by  adjusting  i  and  a  properly,  one  can  obtain  a  good  fit  to  the  data. 

What  factor  was  the  most  important  for  the  success  of  the  model  of  VanZandt 

3 

et  al?  Their  answer  to  this  was  that  it  is  the  variation  of  the  quantity  F.  (the 

fraction  of  turbulence)  which  was  most  responsible  for  the  agreement  between 
2  2 

(C  (model))  and  (C  (radar)).  This  accounted  for  both  the  day-to-day  variation 

arid  the  variation  with  altitude.  Both  types  of  variations  were  found  to  be  very 

large,  that  is,  over  an  order  of  magnitude  in  some  cases, 

2 

The  largest  values  of  C  are  found  at  the  lowest,  altitudes.  To  illustrate  this 

n  32 

we  give  Figure  8  which  is  from  Crane.  As  can  be  seen,  there  is  a  very  steep 

2  afrt 

increase  of  C  as  one  approaches  the  ground.  In  this  respect  it  is  unfortunate 
n  3 

that  the  model  of  VanZandt  et  al  has  not  really  been  tested  at  the  lower  altitudes 

33 

(below  7  km  there  are  problems  obtaining  data  and  below  4  km  none  are  available). 

2 

The  reasons  for  data  failure  where  C  is  measured  at  low  altitudes  are  presented 

n3 

in  the  following  list  by  VanZandt  et  al: 

1.  The  radar  receiver  may  be  overloaded  by  the  much  stronger  echoes  at  the 
lower  heights. 


>!< 

Some  applications  would  not  be  harmed  by  this  however. 

30.  Barat,  J.  (1975)  Etude  experimentale  de  la  structure  du  champ  de  turbulence 

dans  la  mo.yene  stratosphere,  C.  R.  Acad.  Sc.  Paris.  280.  Ser.  B. 

691-693. 

31.  Pond,  S. ,  Stewart,  R.W.,  and  Burling,  R.  W.  (1963)  Turbulence  spectra  in 

the  wind  over  waves,  J,  Atm,  Sci.  20:319-321. 

32.  Crane,  R.K.  (1977)  Stratospheric  Turbulence  Analysis,  AFGL-TR-77-0297, 

AD  A 04 7  740. 

33.  Shannon,  R.R.,  Smith,  W.  Scott,  et  al  (1978)  Atmospheric  Phase  Distortions. 

AFGL-TR-78-0178,  AD  A059  402. 
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Figure  8.  Crane's  Data  on  Cn  via  Radar 


2.  The  effect  of  the  separation  between  radar  site  and  location  of  radiosonde 
should  be  larger  at  lower  heights. 

3.  The  radiosonde  measurement  of  humidity  was  quite  poor  in  comparison  to 
that  of  temperature.  At  lowest  altitudes,  the  role  of  humidity  is  dominant.  There¬ 
fore,  the  data  from  the  radiosondes  would  be  expected  to  become  degraded  at 
lowest  altitudes. 

4.  In  addition  there  is  some  evidence  that  the  value  of  i  may  be,  on  occasion 

o 

radically  different  near  the  ground  than  at  higher  elevations,  and  this  is  not  taken 
into  account  by  the  model. 


fvVA 

f  A.:f 


To  the  above  llat  we  may  also  add  the  possibility  that  at  lower  altitudes,  the 
boundary  layer  effect  may  become  Important  (despite  the  fact  that  the  latter  extends 
upwards  for  only  about  1  km  on  average).  If  this  la  correct,  one  would  have  to 
Incorporate  other  parameters  In  addition  to  R^,  parameters  which  are  used  to 
characterize  boundary  layer  turbulence.  The  model  also  leaves  out  effects  of 
convective  storms. 

In  a  private  communication,  Dr.  VanZRndt  mentioned  that  the  following  im- 

2 

provements  will  be  included  in  the  next  C  model  his  group  will  employ: 

n  3 

"The  two  major  conceptual  defects  in  the  model  of  VanZandt  et  al  are,  first, 
the  neglect  of  fluctuations  of  stability,  and  second,  the  use  of  a  constant  value  of 
the  outer  scale  l Q.  The  method  of  correcting  the  first  defect  is  described  in 
VanZandt  et  al.^4  The  second  defect  must  be  corrected  by  introducing  a  probabil¬ 
ity  distribution  of  layer  thicknesses,  but  determination  of  this  distribution  has 
proved  to  be  difficult." 

The  above  will  greatly  reduce  the  need  to  "fit  the  data"  as  one  must  in  the 
present  model  (for  determining  t  and  a).  Since  this  work  is  at  the  forefront,  we 
await  the  results  of  the  new  model  with  interest. 

4.2  Previou*  Work 

32  2 

Crane  has  made  observations  of  Cn  with  the  Millstone  Radar  facility.  In 

addition,  he  too  had  the  problem  of  estimating  the  fraction  of  turbulence 

which  occurs  below  the  scale  of  the  radiosonde  radar  resolution  (again  of  order 

2 

1  km).  In  addition,  he  has  estimated  C,p  profiles,  t  (dissipation  rate)  and 

(local  eddy  diffusivity)  on  the  basis  of  Cj*  profiles  determined  via  radar.  It  is 

useful  to  compare  some  of  the  things  that  he  did  which  might  have  value  as  parts 
2 

of  future  C  models, 
n 

Crane,  who  is  a  pioneer  in  this  field,  made  his  first  contribution  in  this  area 
with  his  Ph.  D  thesis  Microwave  Propagation  Through  a  Turbulent  Atmosphere. 

In  this  he  used  a  formulation  where  e,  in  contrast  to  Eq.  (3  5)  or  to  Eq.  (47),  is 
not  given  by 

e  --  Kn  (S)2  (109) 


34.  VanZandt,  T.E.,  Gage,  K.S.,  Warnock,  J.M.  (1979)  A  statistical  model  for 

the  probability  of  turbulence  and  the  calculation  of  vertical  profiles  of 
turbulence  parameters,  Proc.  Conf.  Fourth  Symp,  on  Turbulence.  Diffus¬ 
ion,  and  Air  Pollution,  Jan.  15-18,  1979,  Reno,  Nev. ,  pp.  Sz-55. 

35.  Crane,  R.K.  (1969)  Microwave  Propagation  Through  a  Turbulent  Atmosphere, 

Ph.D  thesis,  Dept,  of  Electrical  Engineering,  Worcester  Polytechnic 
Institute. 
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(where  S  is  mean  shear).  The  more  accurate  relation  given  by  Crane  is 


S  -  B 


where  B  is  the  rate  of  buoyancy  dissipation  given  by 

B  *  -  (O'  w'>  —  .  (Ill) 

6o 

<w'  is  the  vertical  velocity  fluctuation) 

This  represents  the  rate  at  which  energy  is  given  to  the  potential  energy  cascade. 
Using  the  definition  of  Kg  given  by  Eq.  (16)  in  a  form  appropriate  for  the  present 
situation,  that  is 

K  =~ie_'..w!,)  (112) 

u  (d?/dz) 
we  obtain 


B  -f  Ke  (ij  =KX 


where  Ng  is  the  buoyancy  frequency  defined  by 


..2  H  g 

B  S  ** 
o 


then  from  Eq.  (110) 


«=KMs2-KeNB*KM<s2>  1-K^n, 

M 


Defining  an  "eddy  Prandtl  number"  by  P£  given  by 


pe  ■  Km/Ks 


we  have 


6  =  Si  S2  (1  -  Pe 1  Ri> 


v  *4 


—  1  36 

At  this  point.  Crane  made  a  model  for  based  on  the  work  of  Ellison,  It 


P  _1  =  (1  +  0,  4  R,)'1  1  -  0.  946 

e  i 


tan'1  (R.) 


In  a  report  to  be  published,  '  he  has  simplified  this  to 


P£«  4R(  .  (119) 

I  do  not  know  If  these  refinements,  that  Is,  that  of  taking  buoyancy  dissipation 

Into  account  or  P  Into  account,  are  of  much  significance;  but,  this  can  only  be 

decided  later  when  more  refined  measurements  and  models  have  been  made. 

32 

The  other  Interesting  point  made  in  the  work  by  Crane  is  that  one  can  esti¬ 
mate  the  fraction  of  altitude  which  is  turbulent  in  the  following  way. 

Let  AV  be  the  difference  in  velocity  across  a  given  height  region  as  measured 
by  means  of  a  radiosonde.  The  resolution  of  such  a  measurement  is  expected  to 
be  sufficiently  bad  as  to  be  far  too  smeared  out  to  allow  one  to  find  layers  where 

o 

R.  <  1/4  (as  was  mentioned  in  the  above  descriptions  of  the  work  of  VanZandt  et  al  ). 

1  3 

To  circumvent  this  difficulty  (that  is,  to  estimate  what  VanZandt  et  al  called  F), 

32 

Crane  proceeded  as  follows;  Let  Av  be  the  velocity  differential  across  a  layer 
of  thickness  below  the  radiosonde  height  resolution,  such  that  the  layer  will  be 
unstable  and  become  turbulent.  Let 


Av  =  AV  . 


(120) 


This  Is  simply  a  bald  assumption.  Next,  define  R.  as  the  gradient  Richardson's 
number  based  on  average  vertical  potential  temperature  gradient  and  mean  vertical 
shear  of  the  horizontal  winds,  thus 

F~  =  Ng/[AV/h]  2  (121) 

where  h  is  the  height  resolution  of  the  radiosonde.  Now,  the  Richardson  number 
across  the  actual  turbulent  layer  thickness,  d  (where  d  cannot  be  measured  directly) 
is  given  by 


36.  Ellison,  T.  H.  (1957)  Turbulent  transmission  of  heat  and  momentum  from  an 

infinite  rough  plane,  J.  Fluid  Mech.  2:456-466. 

37.  Crane,  R.K.  (1980)  Radar  observations  of  turbulence  in  the  lower  strato¬ 

sphere,  Radio  Science,  special  issue  (in  press). 
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(122) 


R.  =  - e-s  *  ‘ 

1  (Av/d)2  T 

2 

In  other  words  we  assume:  (a)  N-.  across  d  and  h  are  the  same;  and  (b)  when 

1  “  2 
Ri  =  T’  t*ie  lftyer  wil^  become  turbulent.  By  eliminating  Ng  from  Eqs.  (121)  and 

(122)  we  obtain 

(ff-T^i)'1  (123) 

and  this  is  Crane's  estimate  of  F,  that  is,  (d/n)  is  the  fraction  of  h  that  he  estimates 
to  be  turbulent. 

3 

How  does  this  compare  to  the  method  of  VanZandt  et  al?  On  one  hand  it  is 

simple.  On  the  other  hand,  it  appears  to  be  rather  arbitrary.  With  regard  to 

"success,"  it  has  not  been  tested  in  the  same  manner  as  has  the  other  model 
2  2 

(that  is,  via  (C„  (model))  vs  (C„  (radar))  comparisons).  Yet,  according  to 
37  n  n  — 

Crane  the  average  value  of  F  (call  it  F)  obtained  in  the  above  manner  from  hi3 
Millstone  data  is  (for  the  stratosphere) 

F  -  0.026  (124) 

CRANE 

3 

which  is  very  close  to  that  found  by  VanZandt  et  al,  that  is, 

F  =  0.  03  (125) 

V.  Z. 


(Incidentally,  F  -  0.  1  for  the  troposphere.) 

V.  Z. 

Thus,  the  recent  work  of  Crane  is  in  average  agreement  with  the  very  different 

3 

approach  of  VanZandt  et  al  in  regard  to  F  calculations.  The  two  approaches  thus 
lend  to  each  other  some  support. 


5.  SUGGESTIONS  FOR  FUTURE  RESEARCH 

5.1  Are  There  Geographic  and/or  Seasonal  Dependencies  of  Cf"? 

3  3 

This  is  a  very  interesting  and,  as  yet,  unanswered  question.  Shannon  et  al 
have  said  "the  probability  of  geographic  variation  can  be  considered.  The  four 


sites  used  "(Florida,  Colorado,  California  and  Arizona)"  were  chosen  for  various 

specific  reasons,  however,  the  statistics  of  the  phase  power  spectrum  values  were 

generally  independent  of  geographic  location.  Differences  noted  were  due  to  the 

diurnal  cycle,  but  not  (sic.)  consistent  dtfferentation  due  to  season  or  locality  was 

noted."  In  other  words  they  did  not  see  a  geographic  or  seasonal  dependence  of 
2 

C  .  It  would  be  premature  (on  the  basis  of  one  effort)  to  declare  that  in  fact  these 

results  are  true  in  general.  We  know  that  under  certain  circumstances,  geographic 

2 

location  doe3  indeed  affect  C^.  For  example,  a  telescope  sited  on  a  remote  moun¬ 
tain  has  far  better  "seeing"  than  one  at  low  altitude  and  "looking"  through  the 
atmosphere  over  a  large  city.  But,  such  effects  may  be  due  to  altitude  (and  the 
man-made  heat  in  cities,  etc.)  and  geographic  location  per  se  may  have  no  effect. 
What  does  affect  C3?  That,  of  course,  is  the  main  unanswered  question.  Tlme- 
of-day  effects  are  accepted  at  present  (as  has  been  already  mentioned  more  than 
once  in  this  report).  The  effect  of  synoptic  weather  conditions  is  not  yet  known. 

5.2  What  are  the  Hieheat  Priority  Gapa  in  the  Estimation  Capability 

via  Radiosonde  Measurements? 

p 

From  the  above  it  is  clear  that  the  greatest  deficiency  of  the  C*  model  is  to 
be  found  at  the  lower  height  regions  (below  6  to  8  km)  and  extending  all  the  way  to 
the  ground.  Improvement  of  tests  of  the  existing  model  of  VanZandt  et  al3  to  in¬ 
clude  the  lower  altitudes  is  of  high  priority  due  to  the  fact  that  the  highest  values  of 
2 

are  also  located  in  this  lowest  height  region.  Of  course,  due  to  weighting  effects, 
higher  altitude  C3  can  exert  a  disproportionate  influence.  As  already  has  been 
mentioned,  some  applications  do  not  involve  lower  altitudes.  Nevertheless,  a 
comprehensive  model  would  be  desirable.  As  has  also  been  mentioned,  however, 
work  is  already  in  progress  to  improve  the  model  itself.  (At  present,  however, 
there  is  no  work  being  done  which  takes  into  account  the  boundary  layer  effects  upon 

5.3  Experimental  Suggestion! 

33 

1.  The  report  by  Shannon  et  al  used  laser  sources  mounted  on  a  U-2  air¬ 
craft  and  measuring  devices  on  the  ground.  In  this  way  they  ascertained  optical 

information  about  the  C3  properties  of  the  atmosphere  between  the  source  and 
n  2 

receiver.  The  integrated  values  of  Cn  could  then  be  compared  with  radar  measure¬ 
ments.  Since  this  report  represented  the  only  comparison  between  radar  and 
optical  measurements  that  I  could  find,  it  would  appear  to  be  an  important  source 
of  information  regarding  future  experiments.  One  of  the  most  important  facts  to 
emerge  from  this  report  was  that,  "It  is  in  fact  a  conclusion  of  this  study  that  no 
further  attempts  should  be  made  to  obtain  atmospheric  data  from  thermosondes." 
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Thus,  unless  there  were  an  attempt  to  advance  the  thermosonde  technology  (ad- 

33 

dressed  to  the  problems  encountered  by  Shannon  et  al  )  one  should  avoid  making 
such  measurements.  On  the  other  hand,  it  may  be  possible  to  overcome  the 
problems  mentioned  provided  the  difficulties  can  be  pinpointed  ar.d  appropriately 
handled  without  undue  cost. 

2.  Simultaneous  radar,  radiosonde  and  smoke  trail  measurements  would  be 

3 

desirable.  Since  the  work  of  VanZandt  et  al  has  relied  completely  upon  "micro- 
shears"  or  high  wind  shears  below  the  resolution  of  the  radiosondes,  it  would  be 
important  to  measure  high  resolution  shears  with  the  aid  of  smoke  trails.  Our 
laboratory  has  developed  a  technique  for  measuring  vertical  wind  shears  of  hori¬ 
zontal  winds  in  the  stratosphere  by  means  of  time  lapse  photography  and  compu¬ 
terized  triangulation  methods  in  conjunction  with  smoke  trails.  The  same  tech¬ 
nology  could  be  used  to  obtain  wind  shears  in  the  troposphere  as  well.  While 
rockets  were  previously  needed  in  order  to  lay  the  trails  in  the  stratosphere, 
techniques  involving  balloon-borne  smok  levices  (subsequently  dropped)  could  be 
used  to  measure  the  actual  microshears.  jVbove  the  boundary  layer  (1-3  km)  and 
in  the  absence  of  weather  fronts,  storms,  etc.,  the  troposphere  behaves  in  a  way 

quite  analogous  to  the  stratosphere.  The  radiosonde  and  radar  observations  could 
2 

be  used  to  calculate  Cn  in  the  manner  described  in  Section  4  (or  by  a  more  ad¬ 
vanced  method  if  available  at  the  time)  and  one  could  compare  measured  micro¬ 
shears  to  the  theoretical  microshears.  One  anticipated  problem  is  that  in  the 
troposphere,  eddy-diffusion  rates  would  be  larger  than  in  the  stratosphere.  This 
may  affect  the  time  history  of  the  smoke  trails  in  an  adverse  way.  F’or  example, 
if  the  eddy  diffusivity  were  very  large,  the  trail  might  spread  too  fast  to  allow  an 

accurate  measurement.  This  would  have  to  be  calculated  in  the  design  of  the 

2 

experiment.  Values  of  C n  would  also  be  available  from  measured  microshears 
and  radiosonde  temperature  profiles. 

3.  Low-altitude  radar  measurements  (below  5  km)  which  cannot  be  made 

with  dipole-array  type  radar  instruments  would  be  available  from  steerable  dish 

radars,  such  as  the  one  located  at  Millstone  (Lincoln  Laboratory).  It  would  be  of 

importance  to  actually  make  such  measurements  in  conjunction  with  radiosonde 

2  2 

measurements  to  see  how  the  (Cn  (radar))  and  (C*  (balloon))  actually  compare. 

Since  specific  humidity  is  of  paramount  importance  at  the  lowest  altitudes,  it 
would  be  necessary  to  pay  special  attention  to  the  precision  of  the  humidity-meas¬ 
uring  instruments. 

2 

4.  It  would  be  desirable  to  have  in  situ  measurements  of  C  as  measured 

-  n 

from  a  balloon-borne  device.  Since,  it  is  not  yet  obvious  that  thermosondes  are 
good  fo"  this  purpose,  it  would  be  very  valuable  if  the  precise  causes  (of  the  in-situ 
measurement  problem)  were  ascertained.  It  would  be  an  important  breakthrough 
if  a  technique  were  devised  so  that  such  measurements  could  be  rendered  possible. 
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A  suggested  cause  of  thermosonde  artifact  ts  that  the  balloon  rises  with  such  high 
velocity  that  the  thermosonde  measures  only  wake  effects.  If  this  is  correct,  then 
it  might  imply  that  the  thermosonde  device  itself  might  be  adequate  but  that  a  mod¬ 
ified  balloon  vehicle  would  be  needed.  In  other  words,  it  may  be  necessary  to 
construct  a  small  balloon  of  very  slow  ascent  (or  descent)  in  order  to  make  in  situ 
C  measurements.  Electronic  problems  also  seem  to  exist  as  well  and  these  too 
would  have  to  receive  attention. 

5.  Research  on  the  nature  of  clear  air  turbulence  at  all  altitudes  is  needed. 

2 

This  would  represent  the  only  hope  for  a  method  to  predict  C  variations.  Along 

**  38 

these  lines,  a  report  on  the  nature  of  CAT  was  put  out  by  this  organization 

which  indicated  that  aircraft  in  situ  measurements  of  supposed  turbulence  may 

often  measure  gravity  waves  instead  of  turbulence.  This  may  imply  that  gravity 

waves,  perhaps  trapped  gravity  waves,  are  the  cause  of  the  turbulence  (and 
2 

associated  C  ).  If  this  were  the  case,  then  predictive  capability  may  be  made 
possible  from  better  information  on  the  source  of  gravity  waves.  Radar  observa¬ 
tions  of  CAT  in  the  troposphere  and  stratosphere  over  long  periods  of  time  and 
with  high  resolution  (preferably  of  order  100  m  or  better)  would  be  able  to  answer 
many  of  the  presently  unanswered  questions  along  the  lines  of:  (a)  Are  there  pre¬ 
ferred  altitude  regions  for  turbulence  to  occur?  (b)  How  long  do  these  last? 

(c)  Does  the  C*  pulsate  as  if  the  turbulence  were  due  to  shears  generated  by  a 
traveling  wave  (or  waves)?  (d)  From  the  specular  reflection  of  stable  layers  can 
one  ascertain  if  there  are  regions  of  gravity  waves  without  turbulence?  (e)  From 
radar  measurements  of  velocity,  can  the  wind  profiles  "trap"  observed  waves? 

Many  other  experiments  along  these  lines  are  suggested  by  the  observations 
3  9 

of  Woods  and  Wiley  in  the  upper  ocean  such  as  "are  there  ensembles  of  turbulent 
layers?"  etc.  Perhaps  high-resolution  radar  measurements  are  the  only  feasible 
ones  at  present  to  investigate  such  questions;  however,  balloon-borne  anomometers 
must  also  be  given  due  consideration. 

6.  CONCLUSIONS 

Optical  turbulence  is  a  crucial  factor  in  lhe  design  of  some  of  the  laser  defense 
systems  now  under  development.  An  improved  data  base  is  therefore  needed  in  this 
connection  as  well  as  a  better  physical  understanding  of  the  nature  of  the  phenomenon 
itself. 

38.  Dewan,  E.M.  ( 1979b)  Stratospheric  wave  sped  ra  resembling  turbulence, 

Science  204:832. 

39.  Woods,  J.D.,  and  Wiley,  R.  L.  (1972)  Billow  turbulence  and  ocean  micro- 

structure,  Deep  Sea  Res.  19:87-121. 


This  subject  is  somewhat  urgent  due  to  the  fact  that  turbulence  compensation 
systems,  if  underdesigned,  will  bo  unreliable,  whereas  If  overdesigned  will  Incur 
very  great  but  unnecessary  financial  cost  as  well  as  delay  In  acquisition. 

It  is  hoped  that  the  preceding  pages  will  acquaint  the  reader  with  some  of  the 
issues  on  this  subject  and  Introduce  him  to  some  of  the  literature. 
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Appendix  A 

Mathematical  Approach  of  Tatanki 

o 

In  the  text,  a  physical  approach  was  given  in  order  to  derive  in  terms  of 
eN  and  e.  A  more  basic  and  exact  approach  is  available,  and  in  this  appendix 
the  latter  is  described.  We  shall  be  following,  in  a  simplified  way,  the  treatment 
given  by  Tatarski.  1  Our  objective  will  be  to  arrive  at  the  following  relation 

Kgigrad  0>2  =  v q  «grad  9')2)  (Al> 

where  0,  as  in  the  text,  corresponds  to  an  additive  scalar,  fl'  is  the  fluctuation 
value,  0  the  average  value,  and  (  )  is  also  a  symbol  for  averaging.  Thus 

0  s  0  +  0'  (A2) 

The  starting  point  for  this  treatment  is  the  molecular  diffusion  equation 

^  +  div  (-Vq  grad  0)  =  0  (A3) 

where  i/g  is  the  molecular  diffusion  coefficient  for  0  and 

5F  K  IF  +  v  •  grad  0  (A4) 
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is  the  mobile  derivative.  Assuming  Incompressibility,  which  in  the  present  case 
Is  permitted,  we  take 


Is  the  turbulent  flux  of  0.  Notice  that  we  switch  from  vector  notation  to  component 
notation  as  convenience  dictates  In  this  treatment.  As  In  the  text  we  define  the 
eddy  diffusion  coefficient  by  0  by 


qT  ■  -  grad  0 


Next,  a  measure  of  fluctuation  Inhomogeneity,  G,  Is  defined  by 


G  ■  I  /  <0'2)  dV 


If  (  )  were  zero  everywhere  inside  of  volume  V,  then  G  =  0.  Subtracting  Eq. 

(A  10)  from  Eq.  (A7)  one  can  obtain  the  equation  for  the  time  rate  of  change  of  0'. 
First  we  use  Eq.  (A2)  and  Eq.  (A8)  in  Eq.  (A7)  before  carrying  out  the  subtraction 
in  order  to  decompose  It  into 

§r~~'  '  lr+l^  [<Vvi><?+e,)  •  ve  ixY'"'~'  ■  vi  ^  ■  <vt e ') +  =  0  <A16) 

and,  using  Eq.  (A  10),  that  is,  subtracting  the  average  equation,  we  obtain 


-+fc-[vie+7ifl,+v; e,-<v! e,>  u6 


At  this  point,  we  use  v^  =  v.  +  v.  in  the  coefficient  of  0'  to  simply  Eq.  (A  17) 


If  -o 


We  may  now  employ  Eq.  (A  18)  to  obtain  the  equation  of  motion  for  the  quantity 
2 

(O'  ).  The  average  mean  square  of  the  fluctuation  is  the  crucial  quantity  in  tur¬ 
bulence  studies.  Multiply  Eq.  (A  18)  by  O’ 

e‘  If +  e’  1^7 s->  ■« 

and  employ  the  following  relations  from  elementary  calculus 


0i  M.'  _  iL  (i.  0i2) 

at  at  1 2  ' 


(A  20) 


(A21) 


*'  k  <vi  9'»  ■  k  <».  ?  s,2» 


*'  k  <vi  *  ■  *'  vi  k 


(A22) 


where  v^  is  assumed  homogeneous  (that  is,  independent  of  x^)  in  the  derivation  of 
the  above.  For  reasons  which  will  be  more  apparent  later,  we  need  also  an 
expression  for  fc;  (  Vq  j  which  is  in  the  form  related  to  a  divergence.  This 


la 


v  [ve  HJ]  =  ^7  [}e  01  HJ-]  - 


(A23) 


To  show  this  to  be  correct  one  must  carry  out  the  divergence  on  the  right  hand 
side  and  regard  vq  as  Independent  of  x..  When  Eq.  (A20)  to  Eq.  (A23)  are  all 
substituted  into  Eq.  (A  19),  we  obtain,  resorting  again  to  vectors, 


(^  0'^)  +  dlv  - vq  0'  grad  0'j  +  vq  (grad  0')^  +  9'  v^  grad  0 

-  0'  grad  <y>  O')  =  0 


(A24) 


Eq.  (A 24)  must  now  be  averaged,  and,  in  doing  so  we  make  use  of  a  relation 

[§5E7  <vi  «’>])=  0  <A25> 

This  follows  from 


(A26) 


in  conjunction  with  (0')  =  0. 

Averaging  Eq.  (A24)  with  the  use  of  Eq.  (A25)  we  obtain  the  relation  for 

2 

(01  )  that  we  sought,  namely 


a(0'2) 

at 


+  div 


Vq  (0'  grad  0')  J 


+  (0i  v')  •  grad  0 


+  VQ  ((grad  0')2)  =  0 


(A27) 
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One  integrates  Eq.  (A27)  over  a  volume  V  and  makes  use  of  the  well-known 
fact  that  the  volume  integral  of  a  divergence  is  equal  to  a  surface  Integral  of  the 
normal  flux.  Assuming  the  turbulence  to  be  homogeneous,  the  net  flux  Integrates 
out  to  zero  over  the  surface,  and  therefore 


9G 


sr 


+ 


/ 


dV  =  0 


where  we  have  used 


<A28) 


<«•'  v.)  =  0  (A29) 

in  the  second  term  (of  (A27))  (that  is,  the  total  vector  v  is  replaced  by  v|). 

Recalling  Eq.  (A  14),  the  definition  of  Kg,  and  the  definition  of  qT  given  by 
Eq.  (A  13)  we  obtain 


§f  =  /  lKe(grad  ^  ~  v6  «grad  S')2)]  dV  =  0  (A30) 

V 

Next  we  assume  stationarity  or 

|f  =  0  (A3 1) 

and  continue  to  assume  that  the  turbulence  in  the  volume  we  are  considering  is 
homogeneous  throughout  that  volume.  This  takes  us  from  Eq.  (A30)  to 

Kfl(grad  0)2  =  i/Q  <(grad  0')2>  (A32) 

which  is  the  objective  we  were  seeking. 


Alternatively, 


the  surface  can  be  located  outside  of  the  turbulent  region. 
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Appendix  B 

Diccunion  of  the  Analytic  Approach  to  Inertial 

Range  Turbulence 


In  the  text,  the  commonly  used  dimensional  argument  was  employed  to  derive 
the  r law  for  the  structure  function.  This  leads,  via  an  approach  which  uses 
the  Fourier  transform,  to  the  k~3/3  law  for  the  inertial  range  power  spectrum. 

The  reader,  however,  may  have  mistakenly  received  the  impression  that  no  direct, 
analytic  approach  exists  which  arrives  at  this  result  starting  from  the  Navier 
Stokes  equations  of  fluid  motion.  He  may  also  not  be  aware  of  the  experimental 
evidence  which  has  conclusively  established  this  law.  The  purpose  of  this  appendix 
is  to  elaborate  these  two  facts. 

With  regard  to  the  experimental  evidence,  the  classic  example  is  to  be  found 
in  Grant,  H.  L. ,  Stewert,  R.W.  and  Moillet,  A.  ( 1962)  Turbulence  spectra  from 
a  tidal  channel,  J.  Fluid  Mech.  12:241.  For  recent  references  with  exhaustive 
references  to  early  work  see  Gibson,  Stegun  and  Williams  (1970),  J,  Fluid  Mech. 
4b  153;  Boston  and  Burling  (1972)  J,  Fluid  Mech.  55:473;  and  Wyngaard  and  Pao 
(1972)  Statistical  Models  of  Turbulence,  Springer  Verlag,  Ed.  by  Rosenblatt  and 
Van  Atta,  pp  384-401.  The  refined  version  of  Kolmogorov's  k  5^3  theory  seems 
to  be  experimentally  well  confirmed  at  this  time. 

Regarding  the  analytic  approach,  the  most  successful  is  that  of  Kraichnan;  and, 
a  book  by  Leslie  (1973)  Developments  in  the  Theory  of  Turbulence,  Clarendon 
Press,  reviews  this  approach  in  detail.  The  most  recent  paper  by  Kraichnan  along 
these  lines  seems  to  be  Kraichnan  (1974)  On  Kolmogorov's  inertial-range  theories, 
J.  Fluid  Mech.  62:305-330.  Reviews  will  be  found  in  (a)  Orszag  (1977)  Lectures 
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on  the  statistical  theory  of  turbulence  in  Fluid  Dynamics  edited  by  Ballan  and 
Peube,  Gordon  and  Breach;  and  (b)  Leslie  (1973)  Review  of  developments  in 
turbulence  theory  in  Rep.  Prog,  Phys.  36:1365-1424. 

The  analytic  approach  leads  to  a  constant  of  proportionality  for  the  law 

of  l.  77  which  compares  favorably  with  the  experimental  value  of  1.  5  in  view  of 
the  scatter  in  the  data. 

For  a  one  page  summary  of  Kralchnan's  remarkable  achievement  the  reader 
may  consult  Hinze  (1975)  Turbulence,  McGraw  Hill,  Inc.,  pp  249-250. 


Appendix  C 

A  Guido  to  Further  Reading 

The  following  bibliography  overlaps  some  of  the  references;  but,  the  present 
list  is  singled  out  for  annotation  because  it  contains  key  reviews  with  many  refer¬ 
ences  and  pieces  of  important  information.  In  other  words,  these  would  be  of  high 
priority  to  those  readers  of  this  tutorial  who  wish  to  read  further  on  this  subject. 
The  list  is  biased  very  strongly  in  the  direction  of  atmospheric  physics  rather 
than  optical  physics,  but  some  of  the  reviews  will  nevertheless  help  to  lead  the 
optical  reader  to  the  literature  he  seeks.  All  of  the  articles  are  widely  available. 

A.  Books  and  Reviews 

1.  "  Vertical  profiles  of  refractivity  turbulence  structure  constant:  Comparison 
of  observations  by  the  Sunset  Radar  witha  new  theoretical  model"  by  T.  E.  VanZandt, 
J.L.  Green,  K.S.  Gage  and  W.  L.  Clark  in  Radio  Science  13:819-829  (1978). 

This  is  the  article  referenced  in  Section  4  as  "VanZandt  et  al"  repeatedly.  In 
my  opinion  it  is  the  most  advanced  model  available  in  the  literature  for  relating 
radiosonde  data  to  radar  data.  The  beauty  of  the  model  is  the  combination  of 
success  and  simplicity. 

2.  Wave  Propagation  in  a  Turbulent  Medium  by  Tatarski  (McGraw-Hill,  1961). 
This  is  'the  classic"  on  atmospheric  turbulence  effects  on  radio  wave,  radar  and 
light  propagation.  It  even  includes  effects  on  acoustic  wave  propagation  (which  is 
useful  for  acoustic  C^  measurements  if  one  were  to  consider  them).  This  reference 
is  required  for  anyone  working  in  this  field. 
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3.  Laser  Beam  Propagation  tn  the  Atmosphere  edited  by  Strohbehn 

(Springer-Verlag,  1978).  This  reference  Includes  review  articles  by  several 

2 

experts.  For  lists  of  references  on  for  both  atmospheric  turbulence  and 
optical  effects,  it  is  probably  the  best  and  the  most  up  to  date.  The  present 
report  references  this  book  very  frequently  (for  example,  the  "blob  lens  approach"). 

4.  Laser  Communications  Systems.  Pratt,  W.K.  (John  Wiley  &  Sons,  1969). 
This  book  is  wnltten  at  a  simple  and  concise  level,  and  it  reviews  mainly  the  opti¬ 
cal  and  electronic  sides  of  the  problem.  Someone  who  wanted  to  know  the  "hard- 
ware  side"  of  Cn  dependent  systems  in  a  short  time  would  do  well  to  look  over  this 
book, 

5.  Wave  Propagation  in  a  Random  Medium  by  Chernov  (1960).  This  is  a 

classic  from  Russia  which  is  often  quoted  in  conjunction  with  Tatarski.  It  seems 

2 

to  go  into  more  detail  in  regard  to  the  electromagnetic  side  of  the  picture. 

6.  The  Effects  of  the  Turbulent  Atmosphere  on  Wave  Propagation  by  Tatarski, 
Israel  Program  for  Scientific  Translations,  Jerusalem,  1971  (available  frem  NTIS 
TT68-50464).  This  is  a  more  comprehensive  and  updated  version  of  No.  2  above. 

7.  "Optical  propagation  through  the  turbulent  atmosphere"  in  Progress  in 
Op*'cs,  Vol.  9,  1971  edited  by  Wolf,  Ch.  III.  Also,  Optical  and  Millimeter 
Line-of-sight  Propagation  Effects  in  the  Turbulent  Atmosphere  in  Boundary-Layer 
Meteorology  4,1973,  397-422. 

Both  of  these  review  articles  are  by  Strohbehn  who  edited  reference  no.  3 
above.  He  seems  to  be  the  foremost  reviewer  in  this  field. 

8.  "Consideration  of  atmospheric  turbulence  in  laser  systems  design"  by 

J.  I.  Davis  in  Applied  Optics,  Vol.  5,  pp  139-146.  This  is  the  main  reference 

2 

used  by  Pratt  (reference  4  above)  in  his  one  chapter  on  the  Cn  effects  of  the 
atmosphere.  The  present  report  relied  on  it  heavily. 

9.  Two  issues  of  Radio  Science  have  appeared  which  contain  a  very  large 

number  of  articles  on  atmospheric  turbulence  and  its  effects  on  radio  wave 
2 

propagation  (Cn).  These  are  Radio  Science  4. 1969,  Vol.  12,  and  Radio  Science 
10,1975,  #1. 

These  are  both  special  issues  which  are  exclusively  devoted  to  the  present 
subject, 

1 0 .  Refraction,  Attenuation  and  Backscattering  of  Electromagnetic  Waves  in 
the  Troposphere:  A  Revision  of  Chapter  9,  Handbook  of  Geophysics  and  Space 
Environments,  by  V.J.  Falcone,  Jr.  and  R,  Dyer. 

This  is  an  in-house  report;  however,  tt.-  "handbook"  is  in  the  open  literature 
and  the  number  for  this  report  is  AFCRL  70-0007  Jan.  1970  for  those  who  have 
access  to  DDC  or  who  would  wish  to  write  to  those  authors.  This  report  is  a 
short  and  authoritative  review. 
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B.  Important  Articles 

11.  "The  birth  of  'CAT1  and  microscale  turbulence"  by  D.  Atlas,  J.  Metcalf, 
J.  Richter  and  E.  Gossard  in  Journal  of  the  Atmospheric  Sciences  27, 1970,  pp 
903-913) 

Tins  is  a  classic  because  it  seems  to  be  the  first  published  observation  of 
CAT  (seen  as  CAT)  by  means  of  radar. 

12.  "Structure  of  the  atmosphere  in  the  vicinity  of  large  amplitude  Kelvin- 
Helmholtz  billows,"  by  K.  Browning  in  J.  Roy.  Met.  Soc,  Quart  97:283-299  (1971). 

This  is  a  classic  article  because  it  is  the  first  to  use  a  combination  of  radar 
and  balloon  measurements  to  observe  Kelvin-IIelmholtz  shear-induced  turbulent 
layers.  Their  radiosonde  work  was  of  such  high  resolution  that  they  could  "watch" 
the  R.  number  go  below  1/4  and  then  see  the  turbulence  occur  subsequently  on  the 
radar  (about  one  half  hour  after  R.  went  below  critical).  On  occasion  R^  went 
back  up  above  1/4  before  the  turbulence  had  a  chance  to  develop.  When  the  mixing 
occurred,  the  subsequent  balloon  observations  of  the  temperature  clearly  showed 
the  effect  of  this  mixing.  In  all,  1  7  mixing  events  were  observed. 

13.  "Billow  turbulence  and  ocean  microstructure"  by  J.  P.  Woods  and  R.  L. 
W'iley  in  Peep  Sea  Research,  1972,  pp  87-121. 

I  regard  this  paper  as  the  description  of  an  analogue  computer  model  of  strato¬ 
spheric  turbulence  and  CAT  in  the  tropoSDhere.  A  careful  reading  of  this  paper 
will  suggest  many  things  to  a  person  inter  ested  in  the  nature  of  atmospheric  tur¬ 
bulence. 

14.  "The  Dynamics  of  the  Upper  Ocean"  by  O.  M.  Phillips,  Cambridge  Univ. 
Press,  1977,  2nd  edition.  Chapter  5  on  internal  waves  is  crucial  to  an  under¬ 
standing  of  certain  key  factors  in  atmospheric  buoyancy  (gravity)  waves  and 
turbulence. 

15.  "Turbulence  in  stably  stratified  fluids:  A  review  of  laboratory  experi¬ 
ments"  by  S.  A.  Thorpe  in  Boundary  Payer  Met.  5  (1973)  (p  94-119).  This  is  an 
excellent  review  of  the  subject, 

lli.  "Instability  and  turbulence  in  a  stratified  fluid  with  shear"  by  Koop  and 
Browand  in  J.  Fluid  Mech.  93  (1979)  pp  135-153,  This  is  the  most  recent  work, 

I  think,  on  laboratory  studies  of  the  shear  instability, 

17.  The  following  reports  that  I  have  published  are  of  interest  in  connection 
with  the  physical  nature  of  clear  air  turbulence. 

Rosenberg,  N.  and  Dewan,  K..  Stratospheric  Turbulence  and  Vertical  Effective 
Diffusion  Coefficients,  A FCR L-TR -75-05 19.  This  paper  has  been  quoted  by 
several  authors  whose  work  is  in  the  radar  detection  of  turbulence  such  as 
VanZandt  et  al,  Balsley  et  al,  Crtne,  and  finally  Evans.  It  provides  a  review 
of  clear  air  turbulence  structure  and  physics,  and  it  contains  useful  references. 


Dewan,  E.  "Stratospheric  wave  spectra  resembling  turbulence,"  Science 
204:832-835  (1979).  This  paper  gives  evidence  and  theory  for  the  possibility  that 
waves  and  turbulence  may  be  parts  of  a  single  cascade  phenomenon. 

Of  some  possible  use  for  background  Information  are  "Mixing  In  Billow 
Turbulence  and  Stratospheric  Eddy  Diffusion.  "  AFGL-TR-79-0091,  "Estimates 
of  Vertical  Eddy  Diffusion  Due  to  Turbulent  Layers  In  the  Stratosphere,  " 
AFGL-TR-79-0042,  and  "Theoretical  Explanation  of  Spectral  Slopes  in  Strato¬ 
spheric  Turbulence  Data  and  Implications  for  Vertical  Transport."  AFGL-TR- 
76-0247. 
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